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This  report  presents  a  comprehensive  review  of  the  state  of  the  art  of 
probabilistic  seismic  hazard  analysis.  This  type  of  analysis  essentially  pro¬ 
vides  Information  regarding  the  likelihood  that  various  levels  of  a  ground 
motion  parameter  will  be  exceeded  in  a  given  period  of  time. 

The  report  begins  with  a  discussion  of  the  various  seismic  parameters  re¬ 
quired  for  the  analysis  and  the  calculation  of  these  parameters  based  on 
seismic  history  data.  Problems  Involved  In  the  estimation  of  these 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  In  this  report  can  be  converted  to 


metric  (SI)  units  as  follows: 

Multiply 

By 

To  Obtain 

feet 

0.3048 

metres 

miles  (U.  S.  statute) 

1.609344 

kilometres 

square  miles  (U.  S.  statute) 

2.589988 

square  kilometres 
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to*  ft, M 'm i rg  >«c,to*  cnrst f t« to*  4  f.r ,*  i f 4 to  O'4  *.r*  i't  c4 

4*4l/?i*  fto  *•; S i **1  f,  1  ^  f.h«  »/tor.to<5  to.' i"*".*  Of  future  S*’i- 
*)'  «/«*!»$,  The  ',f  the  ter*  "tottol'.  4 1  %lr  4tol/»1»*  tn  earthquake  ervq i neer- 

|*'4  i  i  *  e/  *  f.(f  a  to  refer  tr,  ^.<rh  an  analyst;  ha;  r.au;e<j  ;o*ne  controversy  In  recent 
year;,  ;1n>e  ;el;<n|r  rl;k  I;  al;o  r}eftne«l  a;  the  eypect^r)  consequences  of  a 
future  ;ol;m|r  eyerif  t>y  planners  arul  ;rime  earthquake  engineers.  In  order  to 
avnM  an/  amhlgulty  In  the  definition,  various  Investigators  (24,  93)  have  adopted 
the  term  "ha/ar d"  to  refer  to  the  probability  of  occurrence  of  a  seismic  event  and 
have  reserved  the  term  "risk"  to  refer  to  the  consequences  of  a  seismic  event, 
this  heporl  will  adhere  tn  Ihe  following  definitions: 

Seismic  Hazard  Analysis:  the  development  of  a  probability  versus  ground 
mol  loo  Intensity  t urve  for  a  given  site. 
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Seismic  Risk  Analysis:  The  development  of  a  probability  versus  consequences 

curve  for  a  given  site.  The  consequences  of  a  seismic  event  may  include 

fatalities.  Injuries,  economic  losses  and  damages. 

Thus,  seismic  risk  analysis,  as  herein  defined,  constitutes  a  complete 
probabilistic  analysis  of  the  future  occurrence  of  seismic  events,  as  well  as  an 
analysis  of  the  expected  consequences  to  a  constructed  facility  of  such  events. 
Thus,  the  results  obtained  from  seismic  hazard  analysis  are  an  input  factor  in 
the  seismic  risk  analysis. 

This  Report  specifically  reviews  the  state  of  the  art  of  probabilistic 
seismic  hazard  analysis  for  a  single  site. 

Benefits  of  Seismic  Hazard  Analysis 

Prior  to  the  development  of  seismic  hazard  analysis  procedures  in  the  late 
1960's,  estimations  of  the  likelihood  of  occurrence  for  seismic  hazards  of  var¬ 
ious  intensities  were  based  on  experience  and  professional  judgment.  Use  of 
seismic  hazard  analysis  procedures  now  provides  greater  logic  and  consistency  in 
such  estimations  as  well  as  more  reliable  information  regarding  the  occurrence 
of  different  seismic  hazards  over  a  given  period  of  time.  The  benefits  can  be 
summarized  as  follows  (23,  35): 

1)  Allows  for  maximum  utilization  in  a  logical  and  consistent  manner 
of  geologic,  geophysical  and  historic  data,  combined  with  the  pro¬ 
fessional  judgment  of  trained  seismologists  and  geologists. 

2)  Accounts  for  uncertainties  in  both  the  data  used  and  the  models 
developed  for  analysis. 

3)  Can  provide  comparative  influences  of  various  parameters  and  as¬ 
sumptions  which  contribute  to  the  overall  computed  hazard,  thus 
Identifying  the  major  parameters  and/or  sources  of  uncertainty 
which  may  require  special  consideration. 
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4)  Provides  crucial  Information  toward  well-balanced  engineering 
designs,  which  require  a  trade-off  between  costly  designs  for 
greater  resistance  and  more  economical  designs  which  run  higher 
risks  of  economic  loss  (23,  50,  51,  88,  104), 

5)  Permits  assessment  of  the  seismic  hazard  In  terms  which  allow 
comparison  with  other  natural  and  man-made  hazards. 

General  Description  of  Methodology 

The  basic  procedure  for  conducting  a  seismic  hazard  analysis  for  a  single 
site  involves  the  steps  outlined  below. 

1)  Local  Seismicity  -  A  study  is  made  of  the  geologic,  tectonic  and  seis¬ 

mic  history  of  the  surrounding  area,  and  potential  earthquake  source  zones  are 
identified.  For  each  of  the  potential  sources  a  prediction  is  make  regarding 
the  future  activity  within  the  source  in  terms  of  expected  number  of  earthquake 
occurrences  of  different  magnitudes  within  a  given  period  of  time.  At  this 
point  in  the  study,  the  probability  that  an  earthquake  of  a  certain  magnitude 
will  occur  within  the  region  under  study  can  be  calculated  by  employing  a  pro¬ 
bability  model  to  describe  the  random  nature  of  earthquake  occurrences. 

2)  Regional  Seismicity  -  An  evaluation  of  the  regional  seismicity  involves 

the  computation  of  the  probability  distribution  of  earthquake  intensities  or 
other  seismic  parameters,  such  as  peak  acceleration,  spectral  velocity,  or  dur¬ 
ation,  at  the  site  under  investigation.  An  attenuation  law  is  applied  at  this 
point  to  describe  the  earthquake  at  the  site  as  a  function  of  earthquake  magni¬ 
tude  and  distance  from  the  site.  For  each  source  zone,  the  probability  of  a 

future  earthquake  exceeding  a  given  value  of  the  specific  seismic  parameter  be¬ 
ing  used  at  the  site  is  computed,  accounting  for  the  random  nature  of  seismic 
activity  and  the  uncertainties  Involved  in  the  estimation  of  the  parameter. 
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The  total  probability  of  exceedence  within  a  given  period  of  time  is  obtained 
by  sunming  the  various  contributions  to  this  probability  of  the  sources  identi¬ 
fied  in  the  first  step. 

3)  Micro-regionalization  -  The  attenuation  law  used  in  the  preceding  step 
very  often  provides  an  estimate  of  the  earthquake  parameter  based  on  firm  ground. 
The  influence  of  the  local  soil  properties  anJ  the  details  of  the  local  geology 
upon  the  computed  parameters  for  firm  ground  can  be  incorporated  in  this  step 
to  arrive  at  a  seismic  hazard  curve  which  relates  a  suitably  chosen  seismic 
parameter  describing  earthquake  shaking  at  the  site  to  the  probability  of  its  ex¬ 
ceedence  in  a  given  number  of  years. 

Scope  of  the  Report 

This  Report  is  limited  to  the  review  of  the  state  of  the  art  of  seismic 
hazard  analysis  for  a  single  site.  Part  2  of  the  Report  discusses  the  various 
sources  of  information  with  regard  to  local  and  regional  seismicity.  Detailed 
descriptions  are  presented  of  the  different  parameters  employed  in  the  analysis, 
as  well  as  the  uncertainties  and  problems  involved  in  their  estimation.  A  dis¬ 
cussion  regarding  the  application  of  Baye's  theorem  in  order  to  combine  data 
from  various  sources  is  also  included.  Part  3  is  a  presentation  of  the  mathe¬ 
matical  formulations  required  for  the  seismic  hazard  analysis,  including  a 
brief  introduction  to  some  of  the  probability  concepts  used  in  such  an  analysis. 
Part  4  deals  with  the  various  applications  of  a  seismic  hazard  analysis  and 
presents  the  recently  published  seismic  hazard  maps.  Part  5  concludes  the 
Report  with  a  summary. 
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PART  2:  SEISMIC  HAZARD  ANALYSIS— INPUT  PARAMETERS 


The  probabilistic  assessment  of  maximum  earthquake  intensities  which  may 
occur  at  a  given  site  requires  input  from  a  number  of  sources  including  geologic, 
geophysical  and  historic  data  regarding  the  seismicity  of  the  area.  None  of 
these  sources  considered  alone  provides  adequate  information  toward  a  reliable 
probabilistic  assessment  of  seismic  hazard. 

It  is  clear,  for  example,  that  the  use  of  historic  seismicity  data  as  the 
basis  for  prediction  of  future  seismic  activity  in  a  certain  region  is  a  highly 
uncertain  and  unreliable  method  (3,  8,  43).  The  use  of  such  historic  data  is 
problematic  in  that  the  available  information  is  often  quite  scarce  and  the  re¬ 
cords  are  extremely  recent  in  terms  of  geologic  time.  Furthermore,  the  infor¬ 
mation  extracted  from  historic  records  often  seems  to  contradict  geologic  evi¬ 
dence  (43). 

Allen  (1976)  has  detailed  the  difficulties  in  estimating  future  seismicity 
on  the  basis  of  short  historic  records.  Historic  seismic  records  for  Eastern 
Mediterranean  regions  and  China,  which  date  back  between  2,000  to  3,000  years, 
show  periods  of  up  to  several  centuries  of  quiescence  alternating  with  equally 
lengthy  periods  of  intense  earthquake  activity.  It  is  apparent  that  estimates 
of  the  frequency  of  seismic  activity  based  on  short-term  historic  data  relative 
to  the  geological  time  scale  are  highly  uncertain.  Thus,  in  regions  such  as 
the  Eastern  United  States  where  historic  data  is  available  for  less  than  a  300- 
year  period  and  does  not,  therefore,  allow  observation  of  a  cyclical  seismic 
energy  release  as  observed  in  Asia,  geologic  and  geophysical  information  is 
most  valuable  in  making  predictions  of  future  seismicity.  Baye's  theorem,  which 
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will  be  detailed  later  In  this  section  of  the  Report,  provides  a  rational  pro¬ 
cedure  for  combination  of  the  various  Information  sources  to  arrive  at  the  most 
reliable  estimate  of  probability  of  future  earthquake  occurrence. 

The  remainder  of  Part  2  will  Include  a  discussion  of  the  parameters  In¬ 
volved  In  seismic  hazard  analysis  as  well  as  the  estimation  of  these  parameters 
based  on  selsmologlcal  data.  Further  discussion  will  deal  with  relevant  geologic 
and  geophysical  Information,  followed  by  a  brief  review  of  Baye's  theorem  and  an 
explanation  of  Its  application  In  seismic  hazard  analysis. 


Analysis  of  Seismic  History  Data 


A  comprehensive  study  of  the  local  seismicity  of  a  region  Includes  the  Iden¬ 
tification  of  potential  earthquake  source  zones;  an  estimation  of  the  average 
rate  of  earthquake  occurrence  for  each  of  the  Identified  sources;  and  the  estab¬ 
lishment  of  an  earthquake  magnitude-frequency  law. 


Seismic  Sources 

The  Initial  task  of  a  seismic  hazard  study  Is  to  Identify  all  potential 
seismic  sources  In  the  region  of  the  given  site.  Cornell  (1968)  has  outlined 
the  following  three  basic  seismic  source  types,  shown  In  Fig.  2.1; 

1)  Point  Source  -  Characterizes  a  potential  source  of  earthquakes,  close¬ 
ly  concentrated  In  space  relative  to  Its  distance  from  the  site. 

2)  Line  Source  -  Describes  potential  earthquakes  originating  from  a  known 
fault. 

3)  Area  Source  -  Utilized  when  the  occurrence  of  earthquakes  In  the  given 
region  does  not  correlate  with  the  geologic  structure  of  the  area,  or  when  the 
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FIG.  2.1  Definition  of  Earthquake  Sources 
(After  Cornell  1968) 


geologic  structure  Is  unobservable  due  to  deep  overburdens.  In  such  an  area, 
earthquakes  are  assumed  to  occur  randomly  and  with  an  equal  probability  of  occur¬ 
rence  throughout  the  source  area. 

The  Initial  establishment  of  all  potential  seismic  sources.  Including  their 
location  In  relation  to  the  site.  Is  essential  to  the  seismic  hazard  analysis. 

The  use  of  historic  records  can  be  useful  In  this  determination,  provided  that 
the  records  cover  a  sufficiently  large  area.  Geologic  and  geophysical  data  pro¬ 
vide  an  Invaluable  Input  In  this  task,  particularly  In  areas  where  faults  have 
been  the  primary  seismic  source  in  the  past. 

A  typical  source  model  is  shown  in  Fig.  2.2.  Here,  Cornell  and  Merz  (1975) 
have  employed  historic  data  to  establish  eight  area  sources  of  seismic  activity 
in  order  to  predict  the  future  seismicity  of  a  site  In  Boston. 

A  seismic  hazard  analysis  may  contain  more  than  one  source  and/or  source 
type  contributing  to  the  overall  probability  of  a  seismic  event  occurring  In  a 
given  period  of  time.  However,  exceptions  to  the  practice  of  using  multiple 
source  types  are  common  in  the  literature.  For  example,  for  reasons  of  simplic¬ 
ity,  computational  efficiency  and  economy,  Donovan  (1978)  treats  both  area  and 
line  sources  as  sets  of  point  sources,  provided  that  the  dimensions  of  the  area 
or  length  which  the  point  represents  Is  considerably  smaller  than  the  distance 
from  the  equivalent  point  source  to  the  site.  A  typical  source  model  used  by 
Donovan  Is  shown  In  Fig.  2.3,  In  which  the  active  faults  located  In  the  region 
under  study  are  represented  by  point  sources. 

For  similar  reasons,  McGuire  (1976)  treats  faults  as  narrow  area  sources 
of  1  to  2  kilometres  In  width.  Probabilities  calculated  using  this  model  of 
fault  representation,  when  evaluated  with  regard  to  calculations  based  on  more 
accurate  models  which  treat  faults  as  line  sources,  have  indicated  approximately 
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FIG.  2.3  Fault  Map  of  Northern  California  (a)  and 
Computer  Model  of  Regional  and  Fault 
Activity  (b)  Used  In  Risk  Analysis 
(After  Donovan  1978) 


a  9*  error  In  the  estimated  hazard  at  all  levels  of  Intensity  for  a  site  located 
a  distance  of  16  km  from  a  fault  which  Is  100  km  In  length.  Studies  show  that 
this  Inaccuracy  or  error  decreases  as  the  distance  from  the  source  to  the  site 
Increases. 

^yyaswanty  et  al.  (1974)  describe  their  preference  of  using  faults  as  the 
only  seismic  source  to  be  considered,  particularly  In  a  region  such  as  California 
where  active  faults  are  scattered  throughout  the  entire  state.  Lomnitz  (1964), 

In  his  study  of  seismic  risk  In  Chile,  has  employed  a  rather  simple  source  model 
establishing  each  location  of  an  historic  seismic  event  as  a  future  point  source. 
It  Is  necessary  to  point  out  that  the  validity  of  such  a  simple  model  is  highly 
questionable,  particularly  if  the  historic  records  cover  only  a  limited  time 
span.  The  heavy  emphasis  in  this  approach  on  historic  and  instrumental  data  dis¬ 
regards  possible  valuable  input  from  geologic  and  geophysical  information. 

Proper  identification  of  all  potential  earthquake  sources  within  a  region, 
and  most  importantly  those  sources  which  are  located  near  to  a  given  site  of 
interest,  is  an  absolutely  essential  prerequisite  to  the  calculation  of  seismic 
hazard  for  the  given  site.  Cornell  and  Vanmarcke  (1969)  have  shown  that  small 
and/or  frequent  events  originating  from  sources  close  to  the  site  provide  the 
most  significant  information  contributing  to  the  total  computed  hazard  for  the 
site  when  the  hazard  is  described  in  terms  of  peak  ground  acceleration.  Thus, 
the  oversight  of  a  potential  source  close  to  the  site  could  result  in  a  signi¬ 
ficant  understatement  of  the  calculated  hazard. 

Rate  of  Earthquake  Occurrence 

Following  identification  of  potential  seismic  sources,  an  average  annual 
rate  of  occurrence  of  earthquakes,  X(mo),  is  assigned  for  each  source.  Rate  of 
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occurrence  can  be  estimated  by  counting  the  number  of  earthquakes  of  engineer¬ 
ing  Interest  with  a  magnitude  exceeding  an  arbitrary  lower  bound,  mo.  end  divid¬ 
ing  this  by  the  observational  period  In  years  of  the  historic  data.  Very  often 
In  seismic  hazard  studies,  a  lower  bound  magnitude,  nto.  Is  chosen  representing 
a  magnitude  cutoff  point  below  which  earthquakes  are  considered  to  have  neglig¬ 
ible  damaging  effects  and  may  be  disregarded  for  the  purpose  of  the  study. 

Thus,  the  parameter  X(m0)  describes  the  average  number  of  earthquakes  with 
magnitudes  larger  than  mo  occurring  annually  within  a  specific  source.  The  ac¬ 
tual  number  to  occur  In  a  specific  source  over  a  given  period  of  time  Is  random. 

In  order  to  model  this  randomness,  seismic  hazard  analysis  most  conmonly  acts  on 
the  assumption  that  earthquake  occurrence  follows  a  poisson  process,  a  model 
which  assumes  that  earthquakes  are  both  spatially  and  temporally  independent. 

The  general  form  of  the  poisson  distribution  Is  given  by  Eq.  2.1: 

P{n,  X(mo)  )  *  X(m0)t  e‘X*mo*t  n  =  0,  1,  2,  .  .  .  (2.1) 

nl 

where  P(n,  X(mQ)  )  represents  the  probability  of  n  events  with  magnitudes  equal 
to  or  exceeding  ny,  occurring  In  a  time  Interval  of  t;  and  X(m0)  is  the  mean  an¬ 
nual  rate  of  occurrence  of  events  with  m  ny,. 

Use  of  the  poisson  process  to  determine  the  rate  of  earthquake  occurrence  has 
been  shown  to  be  Inaccurate  by  a  number  of  Investigators.  Ferraes  (1967b,  1973) 
tested  the  poisson  model  for  earthquakes  in  Mexico  City  and  concluded  that  the 
resulting  data  Indicated  non-randomness  and  thus  the  historic  earthquakes  In  that 
region  were  actually  not  statistically  independent  events.  Ferraes  (1973)  at¬ 
tributed  this  deviation  from  randomness  to  "the  fact  that  the  energy  dissipated 
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by  most  of  the  large  earthquakes  Is  so  large  that  the  disturbance  by  these  earth¬ 
quakes  affects  the  chances  of  the  occurrence  of  subsequent  discrete  earthquakes." 
Based  on  data  from  southern  California,  Knopoff  (1964)  concluded  that  seismic 
events  In  that  region  registering  a  magnitude  of  3.0  or  greater  did  not  correspond 
to  a  polsson  distribution,  thus  suggesting  a  causal  connection  between  successive 
seismic  events.  Lomnltz  (1966)  has  suggested  that  observed  deviation  from  the 
polsson  process  using  statistics  of  past  earthquakes  may  be  due  to  a  lack  of  spa¬ 
tial  rather  than  temporal  randomness.  Isacks  and  Page  (1968),  In  their  comnents 
regarding  the  data  presented  by  Lomnltz  (1966),  suggest  that  the  average  rate  of 
occurrence  of  seismic  events,  X,  Is  a  random  variable,  and  that  discrepancies  re¬ 
garding  the  polsson  process  may  be  attributed  to  this  variation  of  X  with  time. 

As  previously  stated,  variations  of  X  with  time  are  observed  in  historic  data 
of  Eastern  Mediterranean  countries  and  In  China,  where  seismological  records  span 
2,000  to  3,000  years. 

Esteva  (1975)  describes  the  distribution  of  waiting  times  between  earthquakes, 
using  geophysical  information,  by  a  gamma  function,  of  which  the  polsson  process 
Is  a  specific  case.  A  comparison  of  the  expected  costs  of  failure  for  the  two 
processes  has  shown  significant  differences.  Using  the  poisson  model,  the  ex¬ 
pected  cost  of  failure  Is  greater  than  that  for  the  gamma  function  for  a  short 
time  following  the  occurrence  of  an  earthquake.  As  the  period  of  time  Increases 
with  no  further  occurrence,  the  computed  risk  based  on  the  gamma  distribution 
Increases,  eventually  to  twice  the  risk  computed  for  the  poisson  process. 

Although  It  Is  presently  very  difficult  to  substantiate  any  assumption  that 
the  occurrence  of  seismic  activity  follows  a  poisson  process.  It  is  equally  dif¬ 
ficult  to  dismiss  the  assumption  on  the  basis  of  what  little  statistical  data 
exists.  Thus,  the  polsson  distribution  is  frequently  employed  In  seismic  hazard 
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and  seismic  risk  studies.  Rosenblueth  (1969)  suggests  that  the  polsson  assump¬ 
tion  may  be  an  acceptable  model  In  the  design  of  structures  with  life-spans  of 
several  decades.  However,  the  non-polsson  nature  of  seismic  activity  should  be 
carefully  considered  In  decision-making  inmediately  following  a  destructive 
earthquake. 

Magnitude-Frequency 

Two  separate  steps  are  Involved  In  the  probabilistic  estimation  of  future 
earthquake  occurrence  of  a  specified  magnitude,  m.  The  Initial  step  is  a  cal¬ 
culation  of  the  probability  of  occurrence  of  seismic  activity  within  a  given 
period  of  time,  and  the  second  step  Is  the  estimation  of  the  probability  that 
should  an  earthquake  occur  It  will  reach  or  exceed  a  specified  magnitude,  m. 

The  first  probability  calculation  can  be  made  utilizing  the  polsson  model  and 
the  earthquake  occurrence  rate.  A,  as  previously  discussed.  The  second  proba¬ 
bility  estimation,  commonly  referred  to  as  the  conditional  probability,  requires 
an  earthquake  magnitude  vs  frequency  law,  which  establishes  a  relationship  be¬ 
tween  the  number  of  earthquake  occurrences  and  their  respective  magnitudes.  The 
law  most  commonly  applied  in  seismic  hazard  studies  for  this  purpose  is  that  sug¬ 
gested  by  Gutenberg  and  Richter  (1954),  which  assumes  the  existence  of  an  expo¬ 
nential  relationship  between  the  number  of  earthquakes  and  their  respective  mag¬ 
nitudes,  as  shown  in  Eq.  2.2: 

A(m)  «  10a  e'6"1  (2.2) 

This  recurrence  relationship,  illustrated  on  a  semi-log  plot,  yields  a  straight 
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line  given  by  Eq.  2.3: 


log10A(m)  -  «  -  bm  ;  b  •  B/1nlO  (2.3) 

where  A(m)  refers  to  the  mean  number  of  earthquakes  with  magnitudes  equal  to  or 
exceeding  the  specified  magnitude,  m,  occurring  within  a  unit  time;  a  Is  a  con¬ 
stant  that  Is  dependent  on  the  location  and  time  of  the  earthquake  sample  used; 
and  b  represents  a  constant  thought  to  be  characteristic  of  a  particular  re¬ 
gion  (54). 

Assuming  that  there  is  an  equal  likelihood  of  occurrence  for  future  earth¬ 
quakes  along  a  fault  of  length  1,  or  over  a  source  with  an  area  A,  then  Eq.  2.3 
can  also  be  written  as  follows: 

log10X(m)  -  a  -  b  m  (2.4) 

where  A(m)  gives  the  number  of  earthquakes  with  magnitudes  equal  to  or  greater 
than  m  per  unit  area  or  length,  and  per  unit  time;  J  Is  equal  to  a  -  log 
(A  or  1  x  period  of  record);  and  F  «  b  ■  0/lnlO. 

The  recurrence  relationship  of  Eq.  2.3  or  Eq.  2.4  can  be  used  in  order  t: 
determine  the  conditional  probability  of  occurrence  of  an  earthquake  having  or 
exceeding  a  magnitude  of  m. 

If  A(m0)  designates  the  total  number  of  earthquakes  occurring  within  a  unit 
time  and  exceeding  magnitude  iq,,,  then  the  ratio  A(m)/A(m0)  gives  the  probabil¬ 
ity  of  a  future  earthquake  having  or  exceeding  a  magnitude  of  m. 

For  the  purpose  of  Illustration,  consider  the  example  of  a  region  with  his¬ 
toric  records  documenting  the  occurrence  of  twenty  earthquakes  with  magnitudes 
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equaling  or  exceeding  m0  =  4.0,  and  further  documenting  that  five  of  the  twenty 
earthquakes  either  reached  or  exceeded  a  magnitude  of  6.0  on  the  Richter  scale. 
Based  on  historic  seismicity,  the  probability  that  a  future  earthquake  in  the 
region  will  have  a  magnitude  equaling  or  exceeding  6.0  can  be  expressed  by  the 
ratio  5/20. 

Therefore,  the  complementary  cumulative  distribution  function  for  earth¬ 
quake  magnitude  (^(m)  can  be  expressed  by  Eq.  2.5: 

Gp|(m)  =  P[M  >  m|  E]  =  1  -  FM(m)  =  xjm)  =  X(m) 

X(m0) 

=  e-6(m  -  mo)  (2.5) 

where  P[M  >_m|E]  is  the  conditional  probability  of  exceeding  m  given  that  an 
earthquake  will  occur.  The  probability  density  function  for  M  is  given  by 
Eq.  2.6: 

fM(m)  =  dFM(m)  *  Be‘P(m  '  mo)  (2.6) 

din 

Frequently,  the  seismic  history  data  available  for  a  region  is  given  in 
terms  of  earthquake  intensity.  In  such  cases,  the  establishment  of  a  magnitude- 
frequency  law  can  be  accomplished  by  writing  the  recurrence  relationship  of 
Eq.  2.2  In  terms  of  earthquake  Intensity  as  follows: 

X(I)  =  10V6lI°  (2.7) 
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where  I0  represents  the  epicentral  Intensity;  and  c  and  Bj  are  parameters 
which  are  comparable  to  a  and  B  in  Eq.  2.2.  The  parameter  Bi  can  be  related 
to  B,  provided  that  there  is  an  assumed  relationship  between  m  and  I0.  The 
most  commonly  used  equation  regarding  this  relationship  is  that  given  by  Gutenberg 
and  Richter  (1954); 


m  =  1  +  2/3  I0  (2.8) 

Other  reported  relationships  expressing  m  in  terms  of  I0  are  shown  in 
Fig.  2.4.  Using  Eq.  2.8,  the  constant  Bj  can  be  related  to  8  by: 

B  =  3/2  Bi  (2.9) 

The  value  of  the  parameter  B  for  a  particular  region  can  be  estimated 
based  on  data  giving  magnitudes  (or  intensities)  and  frequencies  of  past  earth¬ 
quakes.  For  each  specified  magnitude  m,  the  number  of  earthquakes  with  magni¬ 
tudes  equal  to  or  exceeding  m  is  counted  and  plotted  as  a  function  of  m. 

When  this  number  is  then  divided  by  the  time  length  of  the  record  and  the  area 
it  covers,  the  recurrence  equation  will  be  given  as  in  Eq.  2.4.  The  value  of  B 
is  always  equal  to  the  slope,  b,  of  the  best  fit  line  in  the  semi-log  plot 
multiplied  by  lnlO.  Fig.  2.5  shows  typical  earthquake  magnitude  distributions 
for  the  United  States. 

Alternatively,  the  seismic  history  data  can  be  presented  by  plotting  the 
ratio  X(m)/X(mQ)  vs  m.  The  resulting  relationship  gives  the  probability  of 
exceedence  of  m  vs  m.  The  slope  of  the  best  fit  line  on  a  semi -log  plot  will 
be  equal  to  8,  as  shown  in  Eq.  2.5. 
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MMI  (EPICENTRAL),  IQ 


CURVE 

EQUATION 

REFERENCE 

(1) 

»  ■  1  ♦  1  i„ 

Gutenberg  &  Richter  (1954) 

(2) 

m  =  1.2  +  0.6  IQ 

Northeastern  U.S. 

Chinnery  &  Rogers  (1973) 

(3) 

m  =  2.1  +  0.5  IQ 

Western  U.S. 

Krinitzsky  &  Chang  (1977) 

(4) 

m  =  0.82  +  0.69  IQ 

Washington  &  Oregon 
Algermissen  (1969) 

FIG.  2.4  Earthquake  Magnitude  and  MMI  Relationships 
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NUMBER  OF  EARTHQUAKES  PER  YEAR  PER  SQUARE  MILE,  *  (m) 


4.0  4.5  5.0  5.5  6.0  6.5  7.0  7.5  8.0  8.5  9.0 

EARTHQUAKE  MAGNITUDE,  m 

FIG.  2.5  Earthquake  Distribution  In  the  United  States 
(After  Hsleh  et  al.  1975) 
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This  procedure  can  be  applied  to  earthquake  intensities  as  opposed  to  mag¬ 
nitudes  where  necessary,  as  In  cases  where  the  available  historic  records  are 
given  in  terms  of  earthquake  intensity.  The  computed  slope  Is  then  given  by  p,, 
which  can  be  related  to  6  by  Eq.  2.9.  Typical  complementary  cumulative  functions 
in  terms  of  Intensity  are  given  in  Fig.  2.6. 

The  parameter  6  or  Pj  has  recieved  considerable  attention  from  seismologists, 
geologists  and  earthquake  engineers.  This  parameter  is  significant  in  that  it 
defines  the  relative  probability  of  earthquakes  of  various  magnitudes  occurring 
in  a  specific  region.  The  p  value  is  in  an  Inverse  relationship  to  the  degree 
of  probability,  so  that  the  likelihood  of  a  larger  magnitude  earthquake  occur¬ 
ring  diminishes  as  the  p  value  increases.  The  value  of  P  varies  from  region  to 
region. 

Fig.  2.7,  which  is  reproduced  after  Newmark  and  Rosenblueth  (1971),  gives 
typical  magnitude  recurrence  relationships  for  three  separate  regions.  The 
three  macrozones  represent  three  areas  of  the  earth (s  crust  characterized  by 
different  levels  of  seismic  activity.  The  vertical  axis  in  Fig.  2.7  gives  the 
number  of  earthquake  occurrences  expected  annually  in  each  macrozone.  Tables 
2.1  and  2.2  give  other  published  p  and  values  for  regions  of  the  United 
States  and  other  parts  of  the  world. 

Some  investigations  (28,  47.  58,  79,  81.  102,  103)  have  tested  the  validity  of 
the  recurrence  relationships  suggested  by  Gutenberg  and  Richter  (1954),  given 
in  Eq.  2.2.  Objections  to  the  simplicity  of  such  a  law  have  caused  a  number 
of  modifications  to  be  Introduced. 

A-  Upper  and  Lower  Bounds 

Typically,  p  values  which  are  based  on  empirical  data  range  from  1  to  3. 

Esteva  (1969)  discusses  the  implications  involved  in  using  p  values  of  less 


24 


TABLE  2.1 

VALUES  OF  THE  MAGNITUDE-FREQUENCY  PARAMETER,  0 


REGION 


£  REFERENCE 


Western  Nevada 
Western  United  States 
California 

San  Jose1,  California 
Eastern  United  States 
Southern  New  England 

New  Jersey 

Central  Mississippi 
River  Valley 
Mexico  City 
Circumpacific  Belt 


Alplde  Belt 
Low-Seismicity  Region 
of  World 
Japan 
New  Guinea 
New  Zealand 
Western  Canada 
Central  America 
Columbia- Peru 
Northern  Chile 
Southern  Chile 
Mediterranean 
Iran-Turkmenia 
Java 

East  Africa 


2.1  Douglas  and  Ryall  (1975) 

1.17  -  2.00  Hsieh  et  al.  (1975) 

2.07  Housner  (1969) 

1.2  -  1.3  Donovan  (1978) 

1.93  Hsieh  et  al.  (1975) 

2.19(*0.12)  Chinnery  and 

Rogers  (1973) 

2.17  Isacks  and 

Oliver  (1964) 


2.00(±0.25) 

1.27 

2.16 


1.70 


Nuttli  (1974) 
Ferraes  (1967a) 
Newmark  and 
Rosenblueth 
(1971) 


2.88 

2.81 

3.1 

2.4 

2.5 
3.34 
2.55 
2.0 
2.12 
2.5 
2.7 
2.16 
2.00 


II  II  II 

Dowrick  (1977) 

II  II 

II  It 


II  II 

II  II 


II  II 

II  II 


II  II 

II  II 


II  II 

II  II 


II 


II 
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TABLE  2.2 


VALUES  OF  THE  INTENSITY-FREQUENCY  PARAMETER, 


REGION 

5l 

REFERENCES 

New  York  State  and 

Surrounding  Area 

1.0 

Decapua  and 

Liu  (1974) 

Northern  United  States 

1.05 

Cornell  and 

Merz  (1975) 

Boston  Area 

1.10 

ii  H 

Southern  New  England 

1.31  (±0.07) 

Chinnery  and 
Rogers  (1973) 

Southeastern  United  States 

1.36 

Bollinger  (1973) 

New  Madrid  Zone 

1.43 

Mann  and 

Howe  (1973) 

Mississippi  Valley  - 
St.  Lawrence 

1.17 

Algermissen 

(1969) 

Mississippi  Valley 

0.93 

McClain  and 
Myers  (1970) 

Montana,  Idaho 

1.29 

Algermissen 

(1969) 

Puget  Sound,  Washington 

1.43 

H  II 

Wyoming,  Colorado,  New  Mexico 

1.57 

II  II 

Nebraska,  Kansas,  Oklahoma 

1.13 

II  II 

Oklahoma,  North  Texas 

1.27 

II  II 

California 

1.24 

II  II 

Central  United  States 

1.15 

Liu  and  Fagel 
(1972) 

World 

1.35 

Cornell  and 

Merz  (1975) 

All  Turkey 

1.14 

Gulkan  and 
Yuceman  (1975) 
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than  3  In  the  recurrence  relationship  shown  In  Eq.  2.2.  This  equation,  when 
used  In  conjunction  with  the  relationship  between  magnitude  and  energy  released 
suggested  by  Gutenberg  and  Richter,  predicts  the  liberation  of  Infinite  amounts 
of  energy  In  earthquake  activity  per  unit  time.  If  the  value  of  8  Is  less  than 
3.  For  this  reason,  Cornell  and  Vanmarcke  (1969)  have  proposed  a  modified 
magnitude-frequency  relationship  which  considers  an  upper  and  a  lower  bound. 
Such  a  truncated  linear  relationship  between  the  log  of  A(m)  and  magnitude  m 
is  obtained  by:  1)  assigning  a  probability  of  one  to  the  chance  that  a  fu¬ 
ture  earthquake,  should  It  occur,  will  exceed  the  minimum  magnitude  of  engin¬ 
eering  interest,  m0;  and  2)  assigning  a  probability  of  zero  to  the  chance 
that  a  future  earthquake  will  exceed  a  maximum  "credible"  magnitude,  mj.  The 
proposed  complementary  cumulative  probability  distribution  function  of  mag¬ 
nitude  is  given  by  Eq.  2.10: 

P[M  >  m|E]  =  1  -  FM(m)  =  1  -  (1  -  e~6(m  '  ml>);  m0  <  m  <  m, 

(2.10) 

where 

Kn,1  -  [1  -  e"p(ml  _  mo)]  _1 

Fig.  2.8  shows  the  Influence  of  the  truncated  form  of  the  recurrence  re¬ 
lationship  upon  the  complementary  cumulative  function  for  magnitude.  The  com¬ 
parison  made  in  Fig.  2.8  shows  that  truncation  has  significance  only  if  the  val¬ 
ues  of  6  and  m^  are  both  small. 
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B.  "Nonlinear"  Law 


The  estimation  of  the  value  of  B  using  Eq.  2.2  and  selsmological  data  is 
not  only  difficult,  but  contains  great  uncertainties.  Numerous  examples  are 
present  in  the  literature  in  which  magnitude-frequency  data,  when  plotted  on 
a  semi-logarithmic  plot,  does  not  define  a  straight-line  relationship  over  all 
magnitude  ranges.  Typical  magnitude  distributions  for  the  United  States  were 
shown  in  Fig.  2.5.  The  data  for  the  Western  United  States  clearly  deviates 
from  a  straight  line  relationship  between  log  T(m)  and  m. 

Lomnitz  (1969)  has  equated  the  problem  of  estimating  the  exponent  6  to 
that  of  estimating  the  mean  of  an  exponential  problem.  Data  showing  that  this 
exponential  distribution  is  skewed  indicates  that  the  least  squares  regression 

on  the  data  will  lead  to  an  over-estimation  of  8,  and  the  resulting  estimate 

of  seismicity  will  yield  expected  magnitude  figures  which  are  consistently  low, 
and  therefore  unsafe.  Deviations  of  seismic  data  from  the  exponential  distri¬ 
bution  at  low  magnitude  ranges  are  commonly  observed  and  can  be  accounted  for 
by  the  limited  availability  of  low  magnitude  earthquake  data  in  the  earlier 
years  of  record-keeping,  as  well  as  by  the  uncertainties  present  in  the  data, 
particularly  when  given  in  terms  of  earthquake  epicentral  intensity,  I0.  The 
value  of  8  will  be  under-estimated  when  based  on  data  from  small  magnitude 
earthquakes,  and  thus,  earthquake  frequency  for  larger  magnitude  ranges  will  be 
over-estimated  (see  Fig.  2.5).  In  order  to  correct  for  this,  Ferraes  (1973) 
recommends  using  a  shorter  time  span  and  restricting  the  calculation  to  magni¬ 
tude  ranges  where  the  data  is  reasonably  complete  and  the  recurrence  relation¬ 
ship  is  approximately  a  straight  line  on  the  semi-log  plot. 

Rosenblueth  (1969)  states  that  the  usefulness  of  the  magnitude-frequency 
calculation  is  dependent  on  the  nature  of  the  structure  being  considered  for 
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design.  He  explains  that  the  frequency  calculation  of  large  magnitude  earth¬ 
quakes  Is  most  important  in  regard  to  structures  with  long  economic  life  spans, 
while  the  frequency  calculation  of  smaller  magnitude  earthquakes  will  affect 
the  design  of  structures  with  very  uncertain  strengths.  He  maintains  that  a  re¬ 
currence  relationship  which  Is  obtained  by  a  good  curve  fitting  over  a  wide 
range  of  magnitudes  may  be  adequate  for  some  types  of  structures,  but  that  its 
use  may  not  be  appropriate  in  the  design  of  all  types  of  structures. 

In  an  attempt  to  resolve  some  of  the  difficulties  encountered  in  estimating 
the  value  of  $,  bilinear  frequency  laws  have  been  employed  by  some  investiga¬ 
tors.  Fig.  2.9  shows  a  recurrence  relationship  used  by  Mortgart  et  al.  (1977) 
in  their  seismic  risk  study  of  Costa  Rica.  Two  regression  lines  were  fitted  to 
the  data,  and  a  cutoff  point  for  magnitude  was  introduced  based  on  geologic 
considerations. 

Ferraes  (1967b)  introduced  a  polynomial  to  describe  the  logarithm  of  the 
magnitude-frequency  relationship,  A(m): 

7  3 

log^Q  A(m)  =  bQ  +  b^m  +  b2tn  +  b3m  +  .  .  .  (2.11) 


In  1973,  Merz  and  Cornell,  using  only  the  second  order  polynomial,  intro¬ 
duced  both  lower  and  upper  bound  magnitudes  into  the  magnitude-frequency  law. 
The  resulting  quadratic  law  is  as  follows: 


log1Q  A(m) 


;  m  <  iHq 
2  2 

(m  -  mo)  +  b2  (m  -  m0  ) 

5  m  >  m^ 


m0  <  m  <  rn'i  (2.12) 


31 


F  OCCURRENCES  OF  M  AON!  TUDE  l_  &RCER  THAN  m.  Alfm 


The  cumulative  function  for  m  will  then  be: 


FM(m)  =  P[M^m|E]  =  K^U-expJ 


-  m0)  +  B2(m2  -  n^2) 


]  (2.13) 


where 


K  mi  =  1  -  exp 


- 


m0)  +  62(111! 2  -  m02) 


■) 


81  =  bi  lnlO 


B2  =  b^  lnlO 


Studies  have  shown  that  this  proposed  quadratic  magnitude-frequency  func¬ 
tion  corresponds  well  to  both  regional  and  global  earthquake  data.  Fig.  2.10 
gives  a  typical  magnitude-frequency  function  used  by  Merz  and  Cornell  (1973). 

A  comparison  of  the  hazard  results  based  on  the  quadratic  relationship  given 
in  Figure  2.10  vs  those  based  on  the  "linear"  relationship  is  shown  in  Fig.  2.11. 
From  this  comparison  it  can  be  concluded  that  a  seismic  hazard  analysis  which  is 
based  on  a  "linear"  recurrence  relationship  obtained  by  emphasizing  the  data  at 
small  intensities  may  result  in  an  over-estimation  of  the  hazard. 

It  must  be  emphasized  that  "bilinear"  or  quadratic  formulations  of  mag¬ 
nitude-frequency  laws  should  not  be  indiscriminately  applied  whenever  the  data 
indicates  a  deviation  from  linearity.  Careful  examination  of  the  historic  data 
and  an  evaluation  of  its  validity  and  adequacy  are  essential  prior  to  the  appli¬ 
cation  of  more  sophisticated  relationships  to  account  for  such  deviations. 


NUMBER  OF  EVENTS  GREATER  THAN  OR  EQUAL  TO  MAGNITUOE  m,  X(m) 


FIG.  2.10  Linear  and  Quadratic  Magnitude-  FIG.  2.11  Comparison  of  Hazard  Results 

Frequency  Relationships  (After  Using  Linear  and  Quadratic  Fre- 

Merz  and  Cornell  1973)  quency  Laws  (After  Merz  and 

Cornell  1973) 


AVERAGE  RETURN  PERIOD  IN  YEARS 


All  of  the  preceding  methods  of  estimating  the  probability  that  an  earth¬ 
quake  of  a  specified  magnitude  will  occur  in  a  given  region  and  within  a  given 
period  of  time  have  involved  two  basic  assumptions: 

1)  That  the  number  of  earthquakes  occurring  per  unit  time  follows  the 
poisson  distribution;  and 

2)  That  the  magnitude-frequency  distribution  is  basically  linear  on 
semi-log  paper. 

An  alternate  procedure  for  estimating  the  magnitude-frequency  distribu¬ 
tion  is  based  on  Gumbel's  theory  of  extreme  values  (55).  Based  on  this  theory, 
the  cumulative  distribution  function  of  the  largest  earthquakes  occurring  an¬ 
nually  in  a  specific  region  and  having  magnitudes  less  than  or  equal  to  m  is 
given  by: 

FM(m)  =  e  ae 


(2.14) 

lnf-lnF^m))  =  In  a  -  6m 

The  parameters  a  and  B  are  estimated  through  the  following  process: 

1)  Select  the  largest  annual  earthquake  magnitudes  over  a  period  of  t 
years  and  arrange  them  in  ascending  order. 

2)  Assign  probability  values  F^(m)  for  each  event  using  Eq.  2.15  and  be¬ 
ginning  with  the  smallest  magnitude: 

Fi  (m)  =  -f4  1-1.2 . t  (2.15) 
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3)  Plot  these  probabilities  vs  m  on  Gumbel  probability  paper,  and  es¬ 
timate  a  and  6  based  on  the  best  fit  straight  line. 

Using  the  theory  of  extreme  values,  studies  of  global  data  and  data  from 
southern  California  to  determine  the  probability  of  occurrence  of  maximum  earth¬ 
quake  magnitudes  have  shown  that  the  calculated  probabilities  plotted  fairly 
well  as  a  straight  line  on  Gumbel1 s  probability  paper  (31,  34,  43,  108).  Such 
an  observation  is  not  surprising  and  proves  to  be  mathematically  correct,  assum¬ 
ing  that  earthquakes  occur  as  independent  events  (poisson  arrivals)  and  that  the 
Gutenberg  and  Richter  magnitude-frequency  law  is  valid.  It  is,  in  fact,  possible 
to  relate  the  constants  a  and  b  in  Eq.  2.3  to  a  and  6  in  Eq.  2.14,  as  fol¬ 
lows:  Given  that  the  number  of  years  covered  by  the  record  is  t,  then  the  aver¬ 
age  annual  rate  of  occurrence  for  earthquakes  with  magnitudes  equal  to  or  exceed¬ 
ing  m  can  be  obtained  from  the  Gutenberg  and  Richter  relationship  of  Eq.  2.3. 

By  assuming  poisson  arrivals  of  events,  the  probability  of  at  least  one 
earthquake  exceeding  magnitude  m  occurring  annually  is  given  by: 

P[M  >  m]  =  1  -  e"X(m) 


thus. 


. /_\  ina  -blnlO  m 

P[M  <  m]  =  e”X(m)  =  e-10  e 


(2.16) 


A  comparison  of  Eq.  2.16  with  Eq.  2..14  yields  the  following  relationship: 


a  =  log  (a) 
b  =  0/1 nl 0 


(2.17) 
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Dick  (1965)  has  made  similar  observations  and  has  concluded  that 

"the  statistical  theory  of  extreme  values  gives  formulae  that 
might  very  well  be  expected  from  seismological  evidence,  prin¬ 
cipally  the  frequency-magnitude  law  and  the  fact  that  the  time 
of  occurrence  of  the  largest  earthquake  in  a  six-monthly  per¬ 
iod  tends  to  be  randomly  distributed  in  time." 

Thus,  the  theory  of  extreme  values  provides  an  alternate  method  of  comput¬ 
ing  the  same  parameters  as  computed  in  the  Gutenberg  and  Richter  relationship. 
Further,  it  allows  for  examination  of  the  effects  of  after-shocks  or  for  their 
elimination  from  the  study,  if  necessary  (35).  Donovan  (1978)  states  that,  if 
the  list  of  earthquake  records  is  complete,  the  actual  differences  between  the 
two  approaches  may  be  small.  For  example.  Fig.  2.12  presents  the  results  of 
the  seismic  recurrence  studies  by  Donovan  for  a  region  in  San  Jose,  California. 
This  figure  shows  a  straight  line  relationship  on  Gumbel  paper  established  be¬ 
tween  the  maximum  annual  magnitude,  m,  and  the  probability  that  that  magnitude 
will  not  be  exceeded,  F(m).  The  results  shown  in  Fig.  2.12  were  converted  by 
Eq.  2.17  and  replotted  in  Fig.  2.13  in  order  to  compare  the  recurrence  relation¬ 
ship  obtained  using  the  theory  of  extreme  values  to  the  Gutenberg  and  Richter 
relationship.  The  thiru  line  of  the  figure  is  obtained  by  using  the  maximum 
number  of  events  occurring  during  a  five-year  period.  The  results  are  quite 
similar. 

Upper  Bound  Magnitude 

The  truncated  form  of  the  recurrence  relationship  discussed  earlier  applies 
an  upper  bound  for  earthquake  magnitude.  Often  referred  to  as  the  "maximum  credi¬ 
ble  earthquake",  this  upper  bound  magnitude  is  defined  as  the  largest  earthquake 
likely  to  occur  along  a  seismically  active  fault  in  a  given  region.  The  return 
period  associated  with  such  an  event  is  generally  very  large  as  compared  to  the 
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EXTREME  VALUE  PROBABILITY,  FM  (m) 


EARTHQUAKE  MAGNITUDE.  " 


Extreme  Value  Derivation  of  Seismic  Recurrence  Infor-  FIG.  2.13  Comparison  of  Derived  Re- 

mation  (After  Donovan  1978)  currence  Relationships  for 

San  Jose,  California,  Us¬ 
ing  Gutenberg  and  Richter 
and  Extreme  Value  Tech¬ 
niques  (After  Donovan  1978) 


life  span  of  a  typical  structure. 

The  existence  of  a  finite  upper  bound  has  been  widely  accepted,  and  In  many 
countries  the  design  level  earthquake  shaking  for  such  structures  as  nuclear 
reactors  and  large  dams  Is  based  on  the  maximum  credible  earthquake  Intensity 
(21,  44,  94). 

The  estimation  of  m^  for  a  given  region  or  for  a  potential  earthquake  source 
should  be  based  on  available  geologic  and  seismological  data.  Housner  (1969) 
states  that  low  seismicity  of  a  particular  region  does  not  discount  the  possibil¬ 
ity  of  that  region  experiencing  a  large  magnitude  earthquake,  and  cites  examples 
in  the  United  States  where  large  earthquakes  occurred  in  regions  of  relatively 
low  seismicity  -  specifically  New  Madrid,  Missouri,  in  1811-12;  and  Charleston, 
South  Carolina,  in  1886.  Further  Chinnery  and  Rogers  (1973),  in  their  statistical 
study  of  earthquakes  in  southern  New  England,  argue  that  existing  records  for  the 
region  are  too  recent  to  assume  that  large  magnitude  events  have  not  occurred  in 
the  past,  and  they  consider  the  possibility  that  large  earthquakes  may  occur  in 
the  region  at  infrequent  intervals. 

It  is  clear  that  the  determination  of  m^  cannot  be  based  on  historic  data 
exclusively,  and  that  geologic  data  and  professional  judgment  are  invaluable 
in  such  an  estimation.  However,  the  estimated  value  of  will  always  contain 
a  certain  degree  of  uncertainty.  Esteva  and  Villaverde  (1973)  studied  the  influ¬ 
ence  of  m^  upon  the  computed  annual  probability  of  exceeding  the  peak  acceleration 
and  have  found  that  errors  in  the  estimated  value  of  m^  had  little  influence 
relative  to  high  probabilities  but  showed  significant  influence  at  low  probabil¬ 
ities.  Yegian  and  Whitman  (1978)  have  concluded  that  the  influence  of  upon  the 
computed  probability  of  liquefaction  is  greater  for  smaller  values  of  B.  Fig.  2.8 
shows  recurrence  relationships  for  various  estimations  of  b  and  m^ .  It  is  clear 
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that  the  influence  of  is  most  significant  in  cases  of  a  small  e  value  and 
when  the  design  level  probability  is  low.  The  procedure  for  estimating  m^  sug¬ 
gested  by  Housner  (based  on  seismological  data)  will  be  outlined  and  discussed 
herein.  Other  procedures  based  on  geologic  data  will  be  presented  subsequently. 


Housner* s  Method 

Housner  (1969)  suggests  two  different  approaches  to  the  estimation  of  the 
upper  bound  magnitude  m^ .  The  first  method  is  based  on  a  seismic  zoning  map 
and  involves  the  assumption  that  the  frequency  distribution  curves  for  regions 
of  varying  seismicity  have  different  slopes.  Table  2.3  presents  the  maximum 
earthquake  magnitudes  corresponding  to  each  zone.  Housner  states  that  "these 
values  are  more  or  less  arbitrary  and  are  used  as  a  matter  of  expediency." 

The  second  method  involves  the  assumption  that  the  probability  distribu¬ 
tion  for  m,  ranging  from  zero  to  infinity,  is  the  same  for  all  regions.  This 
assumption  allows  the  maximum  credible  earthquake  magnitudes  for  different 
regions  to  be  related  through  the  use  of  Eq.  2.18: 


Bo  Mn0) 

m2  =  ml  ‘  62  1n(^TmoT) 


(2.18) 


where  m^  represents  the  upper  bound  magnitude;  6  represents  the  magnitude-fre¬ 
quency  parameter  defined  in  Eq.  2.2;  X  is  the  average  seismicity  expressed  in 
terms  of  the  number  of  earthquakes  occurring  per  unit  time;  and  the  numbers  1 
and  2  refer  to  the  two  regions  under  consideration.  Thus,  given  that  the  seis¬ 
micity  for  one  region  is  reliably  determined,  then  the  upper  bound  magnitude  for 
another  can  be  predicted  by  assigning  the  proper  B  and  x  values  for  the  two  re- 
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TABLE  2.3 


ESTIMATED  MAXIMUM  ZONAL  ACCELERATION  AND  MAGNITUDE 
(After  Housner  1969) 


Zone  3  (near  a  great  fault) 

MAXIMUM  ACCELERATION 

50%  g 

MAGNITUDE 

8.5 

Zone  3  (not  near  a  great  fault) 

33%  g 

7.0 

Zone  2 

16%  g 

5.75 

Zone  1 

8%  g 

4.75 

Zone  0 

4%  g 

4.25 
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gions.  Table  2.4  presents  computed  upper  bound  magnitudes  using  California 
as  the  known  region,  where  m^  ■  8.5  and  8^  »  0.48.  Based  on  results  shown 
In  Table  2.4,  one  can  conclude  that  the  upper  bound  magnitude  for  a  region 
with  8  =  0.48  and  with  an  average  number  of  earthquakes  10  times  smaller  per 
unit  area  and  unit  time  will  be  approximately  7.4. 

In  general,  if  only  statistical  data  is  used,  the  maximum  credible  earth¬ 
quake  magnitude  for  a  region  can  be  estimated  based  on  curves  of  Intensity  and 
return  period  and  the  choice  of  acceptable  failure  probability  for  a  given  time 
interval.  Using  a  different  approach,  m^  can  be  estimated  based  on  the  judg¬ 
mental  assessment  of  a  geologist.  Ideally,  both  sources  of  information  should 
be  rationally  combined  (41),  a  procedure  to  be  discussed  in  a  later  section  of 
this  Report. 

Previous  discussion  has  been  limited  to  the  study  of  local  seismicity  and 
to  the  estimation  of  the  likelihood  of  occurrence  of  earthquakes  having  dif¬ 
ferent  magnitudes  or  intensities.  A  complete  seismic  hazard  analysis  for  a 
site  should  further  provide  an  estimation  of  the  probability  that  the  chosen 
earthquake  intensity  parameter  will  exceed  a  specified  value  at  the  site.  In 
order  to  accomplish  this,  an  attenuation  law  is  applied,  essentially  relating 
earthquake  magnitude  or  intensity  (MMI)  and  the  distance  between  the  site  and 
the  epicenter  or  hypocenter  to  an  appropriate  parameter  which  describes  the  in¬ 
tensity  of  the  earthquake  ground  shaking  at  the  site.  The  use  of  an  attenuation 
law  for  estimating  earthquake  shaking  on  firm  ground  constitutes  a  study  of  the 
regional  seismicity.  A  discussion  will  follow  regarding  the  use  of  attenuation 
laws  for  firm  ground,  and  a  later  discussion  will  deal  with  site  effects. 
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TABLE  2.4 


UPPER  BOUND  MAGNITUDE  AS  DETERMINED  BY  SEISMICITY  RATIO 
(From  Housner  1969) 


^X  (^q)/^ 


mi 


1  8.5  (California) 

2  8.15 

4  7.8 


5 

10 

20 


7.7 


7.4 


7.1 


30 


6.9 


100 


6.3 
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Attenuation  Laws 


Earthquake  engineering  utilizes  various  parameters  in  order  to  describe  the 
seismic  hazard  for  a  particular  site.  In  seismic  hazard  analysis,  the  two  most 
widely  used  parameters  have  been  peak  ground  acceleration  and  site  intensity, 
generally  expressed  in  terms  of  Modified  Mercalli  Intensity  (MMI).  Neither  of 
these  provides  an  adequate  description  of  the  severity  of  earthquake  activity  at 
a  given  site  for  most  structural  and  geotechnical  design  purposes.  Recent  stud¬ 
ies  have  attempted  to  consider  other  parameters  in  seismic  hazard  analysis,  in¬ 
cluding  peak  ground  velocity,  displacement  and  spectral  velocity. 

Choice  of  the  proper  attenuation  law  is  a  crucial  step  in  the  seismic  hazard 
analysis.  Numerous  investigations  have  been  conducted  regarding  the  attenuation 
of  earthquake  motions  with  distance,  using  seismological  data.  Table  2.5,  re¬ 
produced  from  McGuire  (1976),  summarizes  most  of  the  published  attenuation  laws. 
While  a  detailed  explanation  of  these  laws  and  their  limitations  is  beyond  the 
scope  of  this  Report,  the  following  brief  discussion  will  point  out  some  of  the 
factors  influencing  the  attenuation  laws  used  in  seismic  hazard  studies. 

A.  Acceleration  attenuation  law 

The  most  commonly  used  attenuation  law  is  one  which  relates  peak  acceleration, 
a,  to  earthquake  magnitude,  m,  and  epicentral  or  hypocentral  distance,  R.  The 
form  of  this  function  is  given  by  Eq.  2.19: 

a  =  b1  eb2m  R‘b3  (2.19) 

where  b^,  b?  and  b^  represent  empirical  constants.  Studies  have  shown  that  appli- 


TABLE  2.5  (Continued) 


Reference 

Data  Source 

Distance  Parameter 

Dependent  Variable 

Equation 

Standard 

Deviation 

Donovan 

San  Fernando, 

Distance  to  energy 

Peak  ground  accel- 

a  *12.783x106(R+25)'2,77 

Not  reported 

0974) 

rock  sites 

San  Fernando, 

center  R  (km) 

. do . 

eration  ag  (gals) 

- do - 

9 

a  =2.054xl05(R+25)*K83 

Not  reported 

soil  sites 

San  Fernando, 

- do - 

. do - 

9 

a  *5.165x105(R+25)*2‘04 

*  ■°-481 
in  ag 

Not  reported 

all  sites 

- do - 

- do - 

9 

ag*1.84xl04  R*1-42 

Western  North 
America,  Japan 
Papua-New  Guinea 

- do . — 

ag*1080  e°-5M(R+25)*1-32 

*  '°'707 
in  a_ 

9 

Donovan 

Worldwide 

Hypocentral  distance. 

Peak  ground  accel- 

a  *1320  e°-58M(R+25)"1’52 

o1n  a  =0.84 
In  ag 

(1973) 

epi central  distance. 

eration  a„  (gals) 

9 

Ouke  and 

San  Fernando, 

or  distance  to 
fault  R  (km) 

Distance  to  energy 

9 

Peak  ground  accel- 

,  6.69  -0.0097R 

an  ~a~ e 

o  *0.052  g 

others 

all  sites 

center  R  (km) 

eration  a„  (gals) 

9  K 

*9 

(1972) 

- do - 

9 

Spectral  accelera- 

S  ,  5J4  e-0.0068R 

o,  *0.053  g 

Esteva  (1970) 

See  reference 

Hypocentral  distance 

tion  sg  (g) 

Peak  ground  accel- 

9  R 

a  *1230  e0-8M(R+25)‘2 

S9 

°ln  ,  *1'2 

!n  ag 

7  °ln  vg*°-84 
(Sheet  2  of  6) 

- do . 

R  (km) 

eration  ag  (gals) 

Peak  ground  velocity 

9 

v  *15  eM(R+0.17e°-59M)*1' 

Vg  (cm/ sec) 

(Continued) 

9 

TABLE  2.5  (Continued) 


Reference 

Data  Source 

Distance  Parameter 

Esteva  and 
Rosenblueth 
(1964) 

West  Coast  of 
United  States 

Hypocentral  distance 
R  (km) 

- do . 

- do— . 

Esteva  and 
Vlllaverde 
(1974) 

Western  United 
States 

- do . 

Hypocentral  distance 
R  (km) 

- . do - 

- do - 

- do - 

Central  United 
States 

. do - 

Gupta  and 

Nuttll  (1975) 

Epi central  distance 
to  Isoselsmal 
a  (km) 

Gutenberg  and 
Richter  (1956) 

California 

Eplcentral  distance 

Housrter  (1965) 

Western  United 
States,  Central 
and  South 

America 

Distance  to  fault 

Dependent  Variable 

Equation 

Standard 

Deviation 

Peak  ground  accel¬ 
eration  ag  (gals) 

a  *2000  e0-8MR'2 

9 

See  reference 

Peak  qround  velocity 
Vg(cm/sec) 

v  *20  eM  R'1  7 

9 

See  reference 

Peak  ground  accel¬ 
eration  ag  (gals) 

ag=5600  e0,8M(R+40)*2 

°ln  a  =0'64 

9 

Peak  ground  velocity 

Vg  (cm/sec) 

vg=32  eM  (R+25)'1,7 

°ln  vg*0'74 

Maximum  average  spec¬ 
tral  acceleration 
x  (gals) 

X=69600  e°-8M(R+70)'2 

°ln  7  '°'75 

Maximum  average  spec¬ 
tral  velocity 

V*250  eM(R+60)'1,7 

°ln  V  =0'64 

V  (cm/ sec) 

Modified  Mercalli 
Intensity 

1*1+3.7-0.0011  a 
e 

-2.7  log  a 

Not  reported 

where  Ifi  Is  epicentral 

intensity 

Peak  ground  accel¬ 
eration 

Graphical 

Not  reported 

Peak  ground  accel¬ 
eration 

Graphical 

Not  reported 

(Continued) 
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TABLE  2.5  (Continued) 


Reference 

Data  Source 

01  stance  Parameter 

Howell  and 

California 

Epi central  distance 

Schul tz 

coast 

to  isoseismal 

(1975) 

a  (km) 

- do — 

Rocky  Mountains, 

Epi central  distance 

Washington, 

to  isoseismal 

Oregon 

a  (km) 

Central  and 

Epicentral  distance 

Eastern  United 

to  isoseismal 

States,  Canada 

a  (km) 

Dependent  Variable 

Equation 

Standard 

Deviation 

Modified  Mercalli 
Intensity 

I=Ie+0. 874-0. 422  In  a 
-0.0186  a 

Oj=0.64 

Logarithm  of  Modified 
Mercalli  Intensity 

In  I=ln  1+0.16-0.0763 
e 

In  a  -  0.0023  a 

Oj=0.44 

Modified  Mercalli 
Intensity 

1=1+1.802-0.628  In  a 
e 

-0.009  a 

Oj=0.61 

Logarithm  of  Modified 
Mercalli  Intensity 

In  1=1 n  1+0.322-0.1098 
e 

In  a  -  0.0012  a 

Oj=0.47 

Modified  Mercalli 
Intensity 

I=Ie+3. 278-0. 989  In  a 
-0.0029  a 

Oj=0.64 

- do - 


McGuire  (1974)  West  Coast  of 
United  States 

- do - 


- do- 


— do - 


Logarithm  of  Modified 
Mercalli  Intensity 

(Other  forms  of  equa¬ 
tions  examined  and 
reported  also) 

Hypocentral  distance 

R  (km) 

Peak  ground  accelera¬ 
tion  ag  (gals) 

Peak  ground  velocity 
Vg  (cm/ sec) 

Peak  ground  displace¬ 
ment  dg  (cm) 

Spectral  velocity 
(T*l  sec;  6=0.02) 
s  (in/sec) 

(Spectral  attenuation 
equations  given  for 
other  periods  and 
dampings  also) 
(Continued) 


In  1=1 n  Ie+0. 480-0. 139  Oj*0.43 

In  a  -  0.00075  a 

where  I»  is  epi central 
intensity 

ag=472xl0°-28M(R+25)*1,3  o1og  #  =0.222 

vg=5.64xl0°-40M(R+25)‘1-2  olQg  y  =0.273 
dg=0.393xl0°-43M(R+25)*0-88  olQg  d  *0.330 
s=0.428x100,38M(R+25)'0-59  olQg  $=0.274 
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TABLE  2.5  (Continued) 


Reference  Data  Source  Distance  Parameter 

Mickey  (1971)  See  reference  Hypocentral  distance 

R  (km) 

- do -  - do - 

- do -  - do - 


Milne  and 
Davenport 
(1969) 

Neumann  (1954) 
Nuttli  (1973a) 

Nuttli  (1973b) 


Western  United 
States,  Central 
America,  Chile 

Eastern  Canada 


West  Coast  of 
United  States 


Central  United 
States 


Central  United 
States 


Epl central  distance 
a  (km) 

- do - 

Epicentral  distance 
4  (mi) 

Epicentral  distance 

Epicentral  distance 


Dependent  Variable 


Equation 


Peak  particle  accel-  a'S.OAxlO0-74"1'^'1  A 
eration  a  (g) 

Peak  particle  velocity  v=4.06xl00-88m'3R"1 '5 
v  (cm/sec) 

Peak  particle  displace-  d^S-BCxlO1 * lm~5R-1 2 


ment  d  (cm) 


Peak  ground  accelera¬ 
tion  ag  (g) 


where  m  is  body-wave 
magnitude 

..  0.69  e1’64" 


,  . 1 . 1M.  2 

1 .1  e  +4 


Standard 

Deviation 

See  reference 

- do - 

- do - 

Not  reported 


Modified  Mercalli 
Intensity 


Modified  Mercalli 
Intensity 


1=1^-9.66-0.0037  4  Oj=0.53 

+1.38  M+0. 00528  4  M 

where  I,  is  site  inten¬ 
sity  '  of  M=7  event 
at  distance  4 

I=Ie+0. 15-3. 17  log  R  Not  reported 
R>1 .12  miles 


Vertical  particle  Graphical  and  tabular  Not  reported 

acceleration,  velo-  for  various  earth- 

city,  and  displace-  quakes 

ment  at  3  frequencies 
for  Rayleigh  waves 

Sustained  ground  accel-  Graphical  and  tabular  Not  reported 
eration,  velocity,  and 
displacement  at  3  fre¬ 
quencies  for  surface 
waves 


(Continued) 
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TABLE  2.5  (Concluded) 


Reference 

Oata  Source 

Distance  Parameter 

Dependent  Variable 

Equation 

Standard 

Deviation 

Orphal  and 
Laboud 
(1974) 

California 

Hypocentral  distance 

R  (km) 

Peak  ground  accel¬ 
eration  ag  (g) 

ag=0. 066x1 00‘4MR* 1,39 

See  reference 

California  and 
nuclear  explo¬ 
sions 

Peak  ground  velocity 

Vg  (cm/sec) 

Vg*0.726x10°*62MR’1-34 

- do - 

- do - 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o 

■o 

1 

1 

1 

1 

1 

1 

1 

I 

Peak  ground  displace¬ 
ment  dg  (cm) 

dg»0.0471xl0°-57MR*1,18 

- do . 

Rasmussen  and 
others  (1974) 

Puget  Sound, 
Washington 

Epi central  distance 

Modified  Mercalli 
Intensity 

Graphical:  data  and 
limits  given  for  each 
earthquake 

Not  calculated; 
data  shown 

Schnabel  and 
Seed  (1973) 

Western  United 
States 

Oi stance  to  fault 

Peak  ground  accel¬ 
eration 

Graphical 

Not  reported 

Stapp  (1971) 

Puget  Sound, 
Washington 

Hypocentral  distance 
to  isoseismal 

R  (km) 

Modified  Mercalli 
Intensity 

I«Ie-2.017  log(R/h) 
-0.008(R/h) 

where  Ifi  is  epicentral 

Not  reported 

intensity,  h  is  focal 
depth 
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cation  of  this  equation  In  cases  of  small  distances  results  In  ground  accelera¬ 
tion  values  which  are  significantly  larger  than  the  observed  values.  For  this 
reason,  the  term  "R"  In  Eq.  2.19  has  been  replaced  with  "R  ♦  25". 

Fig.  2.14  shows  typical  acceleration  attenuation  relationships  for  eastern 
and  western  regions  of  the  United  States.  These  curves,  which  represent  max¬ 
imum  horizontal  acceleration  In  rock,  were  applied  In  the  development  of  a  new 
earthquake  hazard  map  for  the  United  States  by  Algermissen  and  Perkins  (1976). 

The  dashed  curves,  together  with  the  solid  lines  at  close  distances,  correspond 
to  curves  recommended  by  Schnabel  and  Seed  (1973).  The  solid  curves  given  for 
eastern  regions  were  generated  through  the  Introduction  of  various  modifications 
to  the  curves  suggested  by  Schnabel  and  Seed. 

Problems  of  uncertainties  In  acceleration  attenuation  laws  have  t-'en  studied 
by  various  Investigators.  Merz  and  Cornell  (1973)  have  suggested  that  an  error 
term  be  included  In  the  attenuation  law,  as  follows: 

a  -  b1  eb2m  R_b3  •  e  (2.20) 

where  «■  represents  a  random  error  term  lognormally  distributed  with  median  1 
and  standard  deviation  of  .  Donovan  (1973)  has  confirmed  the  lognormal  dis¬ 
tribution  of  the  data  points  around  the  mean.  Values  of  for  different 

attenuation  laws  are  listed  In  Table  2.5. 

The  parameters  b^ ,  b^  and  b^  have  been  treated  as  constants  In  most  atten¬ 
uation  laws.  Donovan  and  Bornsteln  (1970)  have  argued  that  earthquake  magnitude 
and  distance  are  not  Independently  related  to  acceleration,  and  have  conducted  a 
study  In  which  they  have  treated  b  as  a  parameter  dependent  on  distance.  Results 
based  on  data  for  rock  and  hard  soil  Indicated  that  this  parameter  Is  Indeed  de- 
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FIG.  2. 


20  50  100 

DISTANCE  (km) 


14  Acceleration  Attenuation  Curves.  The  solid  lines  are  curves 
used  for  the  eastern  region.  The  dashed  lines  together  with 
solid  lines  at  close  distances  are  the  attenuation  curves 
used  for  the  western  region  (From  Algermissen  and  Perkins 
1976) 


pendent  on  distance.  The  following  are  the  values  obtained  for  b  from  this  study: 


=  2,154,000  (R)'2-10 

b2  =  0.046  +  0.445  log  (R  )  (2.21) 

b3  =  2.515  -  0.486  log  (R) 

where  R  represents  the  distance  to  the  energy  center  or  the  causative  fault  in 
kilometres. 

Fig.  2.15  gives  the  attenuation  relationship  developed  from  these  values. 

The  standard  deviation  was  found  to  decrease  as  acceleration  increased.  This 
variation  is  shown  in  Fig.  2.16. 

B.  Intensity  Attenuation  Law 

As  stated  earlier,  the  Modified  Mercalli  Intensity  (MMI)  is  the  most  commonly 
used  term  for  the  expression  of  site  intensity.  A  description  of  this  scale  is 
given  in  Fig.  2.17.  A  seismic  hazard  analysis  is  performed  in  terms  of  intensity 
when  the  seismological  data  for  the  region  is  given  in  terms  of  intensity,  in  which 
case  an  intensity  attenuation  law  is  required  (26).  Such  a  situation  applies, 
for  example,  to  New  England,  where  most  available  seismological  data  is  based 
on  subjective  interpretation  of  historic  records  of  earthquake  occurrence  and 
damage.  A  second  situation  in  which  results  of  the  hazard  analysis  should  be 
stated  in  terms  of  intensity  exists  when  the  risk  analysis  (which  is  based  on 
the  results  of  the  hazard  analysis)  employs  earthquake  intensities  in  relation 
to  damages.  The  seismic  design  decision  analysis  developed  by  Whitman  et  al .  (1975) 
employs  a  damage  probability  matrix  which  suirmarizes  the  observed  behavior  of 
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STANOAPO  OCVIATKX  (LOO  NO»MAL> 


FIG.  2.15  Attenuation  Relationship  for  Mean  Risk  Acceleration 
at  Rock  and  Hard  Soil  Site 
(After  Donovan  and  Bornstein  1978) 
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FIG.  2.16  Variation  of  Attenuation  Uncertainty 
with  Increased  Acceleration  Values 
(From  Donovan  and  Bornstein  1978) 
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ABRIDGED 

MODIFIED  MERCALLI  INTENSITY  SCALE 


I  Kol  fell  carept  by  a  very  fen  under  especially 

favourable  circumstances. 

n 

Fell  only  by  »  lev  persons  it  real,  ci/w*ci»lly  on 
upper  door*  of  buddings.  Dcliratcly  suspended 
object*  mi y  a«in|. 

m 

FeH  quiic  a«licnb)y  indoor*.  iipciially  nn  upper 
floor*  of  buildinp,  bu*.  many  people  do  nnl  ret* 
ognise  it  M  an  earthquake  Standing  motor  cart 

may  rock  slightly.  Vibration  like  pa**mgof  truck 
Duration  estimated 

nr 

During  the.  day  fell  indoor*  by  many.  outdoor*  by 
few  At  night  »onn  an  aimed  D'*hea.  window, 
doora  disturbed.  «ili»  make  creaking  aound  Sen- 

ration  hkc  heavy  trv»k  striking  budding  Stand- 
ing  molnr  ear*  rorkrd  nolirmldy 

J z 

Felt  by  nearly  ever) one.  many  awakened  Some 
di«ki.  window^,  etc.,  broken,  a  If*  instances  of 
crarked  planter;  unstable  objects  overturned 

Disturbance  of  tree*,  pole*  and  other  tall  objert* 
sometimes  noticed  Pendulum  clock#  mat  atop 

l 

^  Fell  by  all;  many  (tightened  and  run  out  doora 

Some  heavy  furniture  moved,  a  feu  instances  of 
fallen  planter  or  damaged  chimneys.  Damage 
•light 

Hverylindy  run*  outdoor*  Damage  negligible  in 
building*  of  good  dc*»gn  anil  const mciion .  slight 
to  moderate  in  wethhuilt  ordinary  struct  ores. 

enn*iderr»blc  in  poorlv  built  or  had!*  d»**icu*’d 
•Irueiurc*.  *omr  chinincvs  hrokeu  Nm*red  by 
per*onj  driving  motor  r.-ir» 

Damage  alight  in  •|»erii»dv  designed  strueturr*. 
considerable  in  ordinary  suMtantial  building* 
with  partial  collapse,  great  in  poor).*  built  *trur- 
lurea  Panel  nail*  thrown  out  of  frame  structures 

I  all  of  ehiumev*.  lactor*  **.  ntks,  roluniii*.  mouii- 
nientj.  nail*  lli-ivi  furuiiure  overturned  Sand 
and  Mmd  ejer-tr»|  m  «uiall  aniouiif*  <  lianf»»  m 
well  water  Persons  ilru  mg  motor  ear*  di«tmlwd 

1 1 

Damage  considerable  in  specially  designed 
structure*,  well  designed  frame  strueturrs  thrown 
out  of  plumb;  great  io  substantial  building*,  mth 

partial  col  lap**-  Pudding*  shifted  oft  (i«undiiti-*n« 
Ground  cracked  conspicuously.  1‘ndriginund 
pipe*  broken. 

X 

Some  well-built  wooden  Mmcturrs  de%t roved. 
nio*t  nnivmry  aod  frame  *(nir|iirv*  de«t  roved 
with  foundation*,  ground  li.uJiy  rrarkrd  Had* 

>»cnt  band*hdrs  considerable  from  river  haul- 
and  *1eep  slope*  Shifted  mud  and  mud  Wsii-r 
•plashed  (slop*,  ted)  over  banka 

FIG.  2.17  Modified  Mercalli  Intensity  Scale 
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structures  during  earthquakes  as  a  function  of  earthquake  Intensity.  Thus, 
the  risk  studies  described  by  Whitman  et  al.  use  seismic  hazard  results  given  in 
terms  of  MMI. 

The  general  form  of  the  attenuation  relation  can  be  expressed  by: 

!site  =  C1  +  C2  ^epicenter  or  m)  '  C3  1nA 

where  &  represents  the  hypocentral ,  epicentral  or  isoseismal  distance;  and  , 
and  Cg  are  constants  dependent  on  the  geological  features  of  the  area.  In¬ 
tensity-magnitude-distance  relationships  for  the  Eastern  and  Western  United 
States  were  presented  by  Chang  and  Krinitzsky  (1977).  It  has  been  shown  that 
intensity  attenuates  at  a  much  slower  rate  in  the  eastern  portion  of  the  country 
than  in  the  western  regions.  McGuire  (1977c)  has  proposed  relationships  relat¬ 
ing  site  intensity  to  peak  ground  acceleration,  velocity  and  displacement. 

Tables  2.6  and  2.7  give  these  relationships  for  soft  and  medium  sites. 

The  form  of  the  attenuation  law  shown  in  Eq.  2.22  implies  that  the  pre¬ 
dicted  isoseismals  or  lines  of  equal  site  intensities  are  circular.  Hence,  the 
orientation  of  a  potential  earthquake  source  relative  to  the  site  under  study 
is  immaterial,  and  only  its  distance  from  the  site  is  of  relevance.  However, 
in  many  cases,  true  isoseismals  on  firm  ground  are  roughly  elliptical  in  shape, 
and  the  orientation  of  the  corresponding  axes  are  often  correlated  to  the  local 
or  regional  geologic  trends  (44).  Esteva  (1975)  explains  that,  in  some  regions 
such  as  the  eastern  United  States  and  Mexico,  "The  isoseismals  seem  to  elongate 
systematically  in  a  direction  that  is  a  function  of  the  epicentral  coordinates." 
In  such  cases,  the  assessment  of  the  seismic  hazard  at  a  site  should  consider 
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TABLE  2.6 

REGRESSION  COEFFICIENTS  FOR  SOFT  SITES 
(From  McGuire  1977c) 

In  p  »  Cj  +  Cjn  +  C3  In  A  +  C4  I 


(peak  ground  acceleration)  |  2.01 

:ih/sec* 


(peak  ground  velocity) 
cm/ sec 


1  .81 


-1.51 


-1.11 


-1.58 


-1.47 


-.313 


-  .9 


(peak  ground  displacement)  -2.35 
m  -2.67 


.157 

-.398 


In  p 

Ci 

In  ag 

-.831 

(peak  ground  acceleration) 

1.45 

cm/ sec* 

1.47 

in  vg 

-4.02 

(peak  ground  velocity) 

-3.61 

cm/sec 

-3.61 

In  dg 

-4.68 

1 

t/1  j 

| 

in  | 

(peak  ground  displacement) 

cm 

-4.81 

|  -2.67  |  .863  |  -.398 

TABLE  2.7 

REGRESSION  COEFFICIENTS  FOR  MEDIUM  SITES 
(From  McGuire  1977c) 

In  p  =  Cj  ♦  C2m  ♦  Ci  In  A  +  Ct,  ls 


C 


-.359 


-.884 


-.064 


-.776 


.168 


-.509 


In  p 


.781 


In  p 


.753 


the  coordinates  of  the  site  In  relation  to  the  sources. 

Uncertainties  In  Intensity  attenuation  laws  are  caused  by  errors  In  data 
(which  Is  subjectively  recorded),  as  well  as  by  possible  errors  In  the  evalu¬ 
ation  of  the  mean  curve.  Donovan  (1978)  discusses  the  bias  Involved  In  the 
mean  Intensity  attenuation  curve  based  on  Isoselsmal  maps  of  a  region.  A  com¬ 
parison  between  the  mean  attenuation  line  based  on  recorded  Intensities  for  In¬ 
dividual  sites  during  an  event  In  Washington  State  and  the  mean  attenuation  line 
obtained  from  the  Isoselsmals  for  the  same  earthquake  showed  a  difference  be¬ 
tween  the  curves  of  0.73  units  of  MMI  scale,  representing  only  the  probable  error 
In  the  mean  and  not  considering  the  scatter  of  data  about  the  mean.  Cornell  and  Merz 
(1975)  have  used  a  standard  deviation  of  Oj  =  0.5  MMI  for  the  scatter  in  the  mean, 
considering  all  site  conditions.  Table  2.5  gives  a  list  of  standard  deviation 
values  for  various  attenuation  laws. 

C.  Other  Attenuation  Laws  Used 

The  use  of  peak  acceleration  to  describe  the  seismic  hazard  at  a  site  with¬ 
out  due  consideration  to  the  frequency  characteristic  of  the  design  motions  is 
sometimes  Inadequate  for  design  purposes.  Seismic  design  criteria  usually  em¬ 
ploy  a  level  of  design  ground  acceleration  in  order  to  scale  a  set  of  standard¬ 
ized  response  spectra.  This  practice  may  be  adequate  for  the  design  of  short 
period  structures,  but  may  require  modification  in  order  to  accommodate  the  design 
requirements  of  Intermediate  and  long  period  structures.  Cornell  and  Vanmarcke 
(1975)  have  concluded  that  the  response  of  long  period  structures  such  as  off¬ 
shore  drilling  platforms  is  poorly  correlated  to  peak  acceleration  and  suggest 
that  peak  ground  displacement  may  be  a  more  appropriate  parameter  for  use  in 
such  designs.  Therefore,  a  peak  displacement  attenuation  law  may  be  more  useful 
In  seismic  hazard  studies  concerned  with  long  period  structures.  Cornell  and 
Vanmarcke  warn  that  the  use  of  a  peak  displacement  attenuation  law  results  in 
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seismic  hazard  estimates  which  are  more  sensitive  to  the  choice  of  the  upper 
bound  magnitude  than  estimates  which  are  based  on  peak  acceleration.  This  dif¬ 
ference  in  sensitivity  is  attributed  to  the  fact  that  peak  ground  displacement 
is  dependent  on  magnitude  to  a  greater  extent  than  peak  acceleration,  as  Is 
Indicated  by  the  value  of  b^  In  the  attenuation  law.  Further,  while  seismic 
hazard  probabilities  based  on  an  acceleration  attenuation  law  are  sensitive  to 
close  potential  earthquake  sources,  probabilities  based  on  a  displacement  at¬ 
tenuation  law  are  more  sensitive  to  distant  sources  (28). 

Similar  studies  by  McGuire  (1977b)  have  shown  that  the  scaling  of  standard¬ 
ized  response  spectra  for  design  purposes  using  peak  acceleration  together  with 
velocity  leads  to  inconsistencies  regarding  the  calculated  risk  for  medium  to 
long  period  structures.  McGuire  recommends  the  use  of  an  attenuation  law  which 
relates  spectral  velocity,  Sy,  to  earthquake  magnitude,  m,  and  hypocentral  dis¬ 
tance,  R,  as  follows: 

Sy(T)  =  a 1 0bm( R  +  25)  ~c  (2.23) 

The  coefficients  a,  b  and  c  are  obtained  for  each  period  of  a  one-degree-of- free¬ 
dom  system  and  for  a  specific  damping  value  by  applying  regression  analysis  on 
computed  response  spectra  from  recorded  motions.  Table  2.8  gives  the  values  of 
the  coefficients  for  various  frequencies  and  for  2X  damping.  McGuire  (1977a)  has 
also  provided  an  attenuation  function  relating  Fourier  amplitude  spectra  of  hori¬ 
zontal  ground  acceleration  to  earthquake  magnitude,  m,  and  source  to  site  dis¬ 
tance,  R,  as  follows: 
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TABLE  2.8 

SUMMARY  OF  REGRESSION  COEFFICIENTS  FOR  MAXIMUM  PSEUDO- VELOCITY  (IN  INCHES/SEC) 

OF  OSCILLATORS  WITH  Zt  DAMPING 
(From  McGuire  1977b) 


Sv  =  alO^R  +  25)~c 


log'Sv  =  a  4  bM  -  c  log  (R  4  25) 


Frequency,  cps 

a 

t 

a 

O' 

-c 

10 

12.27 

1.051 

.  225 

-1.324 

6.67 

56.89 

1.719 

.185 

-1 . 368 

5 

32.16 

1.476 

.232 

-1.292 

3.33 

29.88 

1.447 

.  295 

-1 . 386 

2.5 

11.42 

1.025 

.321 

-1 . 205 

2 

5.44 

.708 

.360 

-1.146 

1.67 

2.48 

.364 

.440 

-1.203 

1.25 

1.158 

.0318 

.406 

-  .913 

1 

.428 

-  .416 

.382 

-  .587 

0.67 

.1085 

-1.020 

.428 

-  .434 

0.5 

.0918 

-1.100 

.456 

-  .519 

0.33 

.0900 

-1.114 

.473 

-  .615 

0.25 

.0552 

-1.344 

.539 

-  .742 

0.2 

.0735 

-1.223 

.565 

-  .931 

0.167 

.0995 

-1.087 

.465 

-  .694 

0.125 


1825 


802 


413 


705 


(2.24) 


FS(T)  *=  exp  (b,  +  b2m  +  b4Ys)Rb3 

where  FS(T)  Is  the  Fourier  spectrum  ordinate  for  period  T;  Y$  is  a  soil  ampli¬ 
fication  term  and  is  equal  to  one  for  soil  sites  and  zero  for  rock  sites.  The 
coefficients  b^ ,  b2,  b3  and  b^  and  the  standard  deviation  of  InFS  are  giyen  in 
Table  2.9  for  different  values  of  the  period  T. 

Eqs.  2.23  and  2.24  can  be  applied  in  seismic  hazard  analysis  to  compute 
the  annual  probability  of  exceeding  of  Sv  or  FS  vs  the  parameter  Itself. 

Another  parameter  of  significance  in  the  determination  of  earthquake  in¬ 
tensity  or  in  the  specification  of  design  criteria  is  duration  of  earthquake 
motion.  In  geotechnical  problems,  peak  ground  acceleration  and  duration  of 
motion  are  very  often  adopted  to  describe  the  intensity  measure  at  a  site. 

For  example,  it  is  known  that  the  susceptibility  of  a  site  to  liquefaction 
depends  on  the  level  of  acceleration  and  the  duration  of  the  sustained  motion. 

Reported  examples  of  seismic  hazard  analyses  which  have  explicitly  taken 
into  consideration  the  effects  of  duration  are  very  few.  Douglas  and  Ryall  (1975) 
have  described  a  seismic  hazard  analysis  procedure  which  allows  for  the  determin¬ 
ation  of  return  period  of  acceleration  as  a  function  of  magnitude.  Since  dur¬ 
ation  of  earthquake  motion  can  be  correlated  to  magnitude,  the  design  criteria 
corresponding  to  an  acceptable  risk  level  can  be  expressed  in  terms  of  peak 
ground  acceleration  and  duration,  as  in  Fig.  2.18.  A  procedure  which  is  able 
to  explicitly  consider  duration  in  seismic  hazard  analysis  will  require  the 
use  of  an  attenuation  law  for  duration.  The  literature  offers  various  defini¬ 
tions  of  and  attenuation  laws  for  duration  (14,  19,  101).  A  state-of-the-art 
review  of  duration,  spectral  content  and  predominant  period  of  strong  motion 
earthquake  records  has  been  compiled  by  Chang  and  Krinitzsky  (1977).  The  most 
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TABLE  2.9 

SUMMARY  OF  REGRESSION  COEFFICIENTS  FOR  FOURIER  AMPLITUDE  SPECTRUM  FOR  2*  DAMPING 

(From  McGuire  1977a) 


T 

*>1 

*2 

*3 

3>t. 

°Ws 

0.04 

-3.  IB 

0.90 

-3 . 20 

-0.43 

3  .26 

.05 

-3.59 

.75 

-1 .23 

-  .31 

3.37 

.  06 

-0.47 

.73 

-3  -  23 

-  .«£ 

.  95 

.oe 

-0.47 

.9? 

-3  33 

-  .43 

3  .07 

.3 

-0.02 

3 .04 

-3.37 

-  .70 

.  9E 

.13 

-0.62 

.83 

-1.07 

-  .49 

.96 

.17 

3 .06 

.92 

-1.23 

-  .47 

.89 

.2 

-  .35 

1.01 

-1.11 

-  .45 

.92 

.  24 

-  .80 

.99 

-1.19 

-  .12 

.66 

.3 

-  .82 

1.02 

-1.22 

.02 

.86 

.34 

-  .64 

3.05 

-1.36 

-  .07 

.93 

.4 

1.04 

.68 

-  .97 

-  .13 

.  90 

.5 

-  .05 

1.09 

-1.11 

.25 

.93 

.6 

-3.03 

3.16 

-  .96 

.39 

.87 

.8 

-3.03 

3.04 

-  .75 

.36 

.  92 

1 

-1.24 

1.12 

-  .89 

.68 

.93 

3  .3 

-1.84 

1.08 

-  .69 

.60 

.93 

1  .7 

-2.55 

1.32 

-  -85 

.31 

.98 

2 

-3.40 

1.44 

-  .89 

.53 

o 

o 

»-1 

2.4 

-3.08 

1.23 

-  .65 

.  38 

1.05 

3 

-2.54 

o 

o 

n 

-  .46 

.35 

1.03 

3.4 

-2.57 

1.06 

-  .59 

.31 

1.11 

4 

-3.93 

00 

rx 

-  .59 

.30 

1.00 

5 

-3.36 

1.20 

-  .61 

.32 

1  .11 
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DURATION,  D,  IN  SECONDS 


FIG.  2.18  Return  Period  T,  in  years  for  Selected  Values 
of  Rock  Acceleration  a  (in  g's),  as  a  Function 
of  Magnitude  m  and  Duration  D  (in  seconds)  of 
Strong  Shaking  (After  Douglas  and  Ryall  1975) 
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common  definition  of  duration  is  the  "bracketed"  duration  defined  as  the  time 
interval  between  the  first  and  last  peaks  of  acceleration  equal  to  or  greater 
than  0.05  g  for  the  strong  earthquake  record.  Typical  duration  attenuation  laws 
for  soil  and  rock  sites  were  presented  by  Chang  and  Krinitzsky  (1977). 

In  Table  2.10,  duration  is  related  to  site  intensity  for  near  and  far  field 
conditions  and  for  rock  and  soil  sites. 

Mortgart  et  al.  (1977)  have  employed  a  duration  attenuation  law  together  with  an 
acceleration  attenuation  law  in  order  to  describe  the  seismicity  of  a  site.  The 
attenuation  law  used  in  the  study  was  based  on  data  obtained  from  Bolt  (1973). 

The  best  fit  relationship  for  duration  was  found  to  be: 


The  duration  at  50  km  for  various  magnitudes  is  given  in  Table  2.11.  These 
values  were  used  in  Eq.  2.25  to  define  a  complete  duration  attenuation  law  in 
terms  of  m  and  R.  Variance  studies  indicate  that  the  coefficient  of  variation 
was  approximately  0.15. 
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TABLE  2.10 


MM 


DURATION  AS  A  FUNCTION  OF  SITE  INTENSITY  IN  THE  NEAR  AND  FAR  FIELDS 
(From  Chang  and  Krinitzsky  1977) 


NEAR  FIELD  DURATION 

seconds 

FAR  FIELD 

seconds 

DURATION 

OCAL 

1  .fENSITv 

ROCK  SITE 

SOIL  SITE 

ROCK 

SITE 

SOIL  SITE 

V 

6.4 

17.5 

9  to 

10 

17  to 

19 

VI 

7.3 

20.0 

II 

II 

II 

It 

VII 

8.3 

23.0 

II 

II 

II 

II 

VIII 

9.6 

26.5 

II 

II 

II 

II 

IX 

11.0 

30.0 

II 

II 

II 

It 

X 

12.5 

35.0 

II 

II 

II 

II 

XI 

14.5 

40.0 

II 

II 

II 

II 

XII 

16.5 

46.0 

II 

II 

l» 

II 

TABLE  2.11 

BRACKETED  DURATION  (SEC)  AT  50.  KM 

FROM  SOURCE 

(Accel , 

.  -  0.05  g's;  Freq. 

>  2  Hz) 

(From  Mortgart  et  al.  1977) 


Duration 

seconds 

Magnitude 

m 

Duration 

seconds 

Magnitude 

m 

1  -5 

-  0 

17-0 

6.75 

1-^ 

L.?5 

r?.o 

7-0 

1-5 

14-50 

F--?5 

7-?5 

1-5 

^ -75 

?6.o 

7-50 

?.0 

5-0 

?7-?5 

7-75 

?•  0 

5-?5 

OR.O 

6.0 

?•  o 

550 

?8-5 

8.?5 

?.5 

5-75 

?9-o 

6.50 

?.o 

6.0 

?9-  5 

8-75 

5-5 

6.05 

09.5 

9-0 

10.0 

(•50 
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Other  Factors  Influencing  Attenuation  Laws 


A.  Site  Condition 

Incorporation  of  the  Influence  of  the  soil  at  a  site  Is  a  part  of  micro¬ 
regionalization.  Empirical  and  theoretical  procedures  have  been  devised  to 
account  for  this  influence  (4,  5,  12,  26,  39,  46,  77,  90,  91,  92).  In  his  sta¬ 
tistical  analysis  of  the  data  on  the  1971  San  Fernando  earthquake,  Donovan 
(1973)  evaluated  the  differences  between  the  attenuation  characteristics  of 
soil  and  rock.  His  results  indicated  that  the  peak  ground  accelerations  at 
short  distances  on  rock  were  higher  than  those  on  soil.  Results  further  showed 
a  reversal  of  this  trend  at  distances  of  more  than  60  km,  or  at  acceleration 
levels  of  less  than  8%  gravity.  The  resulting  attenuation  law  suggested  for 
soil  and  rock  is  given  in  Table  2.5. 

Another  method  of  incorporating  the  influence  of  soil  conditions  utilizes 
the  curves  by  Seed  et  al,  (1976),  which  are  reproduced  in  Fig.  2.19.  Additional  pro' 
cedures,  which  relate  peak  ground  surface  acceleration  on  soil  to  bedrock  for 
generalized  site  conditions,  are  present  in  the  literature  (8,  17). 

For  a  site  where  there  is  adequate  information  regarding  the  site  and  soil 
conditions,  a  theoretical  analysis  based  on  the  wave  propagation  theory  can  be 
made  to  assess  the  modification  of  the  ground  motions  propagating  from  bedrock 
to  the  surface  of  the  ground.  Faccioli  (1977)  has  described  a  stochastic  method 
of  one-dimensional  soil  amplification  analysis  which  employs  nonlinear  soil 
properties.  Using  this  procedure,  probabilistic  assessments  of  earthquake 
parameters  on  local  soil  sediments  can  be  made  from  the  distributions  of  the 
parameters  for  rock  or  firm  ground. 

Modifications  in  the  common  form  of  the  attenuation  law  have  been  made 
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MAXIMUM  GROUND  SURFACE  ACCELERATION,  g‘s 


0  0.1  0.7  0.3  0.4  0.5  0.6  0.7 


MAXIMUM  ACCELERATION  IN  ROCK,  g't 

FIG.  2.19  Approximate  Relationships  between  Maximum  Accelerations 
on  Rock  and  Other  Local  Site  Conditions 
(From  Seed  et  al.  1976) 
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Introducing  a  suitable  parameter  to  the  equation  In  order  to  represent  the  In¬ 
fluence  of  site  properties.  Performing  multiple  regression  on  the  1971  San 
Fernando  data,  Eguchi  (1977)  has  proposed  Eq.  2.26  for  peak  particle  velocity,  V 


V(cm/sec) 


e4.69k0.33e1.2(m  -  6.4)R-1.05 


(2.26) 


where  R  represents  the  hypocentral  distance  in  kilometers  and  K  Is  a  constant 
which  defines  the  variation  of  shear  velocity  with  depth  of  profile  and  is 
given  by  Eq.  2.27: 

K  =  d0.37  ^2-27) 

where  represents  the  shear  wave  velocity  and  d  represents  the  depth. 

More  recently,  McGuire  (1978)  has  presented  a  modified  form  of  the  atten¬ 
uation  function  in  order  to  incorporate  soil  effects,  as  follows: 

In  p  =  b-j  +  b^m  +  b^  InR  +  b^  Ys  (2.28) 

where  p  is  the  ground  motion  variable  or  the  spectral  velocity;  b^ ,  b^,  b3 
and  b^  are  coefficients  determined  by  regression  analysis;  and  Y$  is  a  site 
geology  indicator  which  Is  equal  to  zero  for  rock  sites  and  one  for  so  sites. 
Table  2.12  provides  the  regression  coefficients  for  peak  ground  acceleration, 
velocity,  displacement  and  spectral  velocity  at  a  period  of  1  second  with  2% 
damping. 

A  similar  study  by  Blume  (1977)  has  resulted  In  site-acceleration-mag- 
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TABLE  2.12 


COEFFICIENTS  AND  STANDARD  DEVIATIONS  OBTAINED  FROM  REGRESSION  ANALYSIS 

(After  McGuire  1978) 


In  p  =  t>!  +  b2  M  +  b3  In  R  +  b4  Yg 


p 

bl 

b2 

b3 

b4 

cln  p 

ag,  peak  ground  acceleration 

3.40 

• 

00 

-1.17 

-  .20 

.62 

Vg,  peak  ground  velocity 

-1.00 

1.07 

-  .96 

+  .07 

.64 

dg,  peak  ground  displacement 

-2.72 

1.00 

-  .63 

.12 

.69 

PSRV,  response  spectrum  ve¬ 
locity  (l-sec  period,  2 
percent  damping) 

-1.61 

1.16 

-  .83 

.31 

.72 

2 

is  cm/sec 
vg  is  cm/ sec 
dg  Is  cm 
PSRV  is  cm/sec 
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nltude  (SAM)  relationships  which  incorporate  the  site  effects.  The  proposed 
expressions  for  the  acceleration  at  any  probability  level,  y,  in  lognormal 
form  are,  for  m  <  6.5: 

Ay  *  0.318e1'03n’(29)1,14E(R  +  25)"1  *14F(2.54)y 

and  for  m  >  6.5:  (2.29) 

Ay  -  26.0e0-432”(29)'22E(R*25)-,-2ZEO.81)i' 

where 

b  =  -y-  log.jQ(pVs) 

in  which  Ay  is  the  estimated  peak  acceleration  that  would  be  recorded  by  an 
instrument  during  an  earthquake  of  m  and  R;  and  pVs  gives  the  site  impedence 
in  feet  per  second.  The  parameter  y  is  the  standardized  variable  of  the  log¬ 
normal  distribution:  when  y  =  0,  the  estimated  acceleration  corresponds  to  the 
mean  value. 

The  Influence  of  soil  conditions  upon  the  intensity  attenuation  laws  is 
more  difficult  to  assess.  The  available  data  on  intensities  corresponds  to  a 
number  of  different  earthquake  locations  and  site  conditions.  This,  coupled 
with  the  fact  that  the  Intensity  scales  are  based  on  subjective  judgments,  ren¬ 
ders  site  influence  determination  nearly  impossible.  Because  intensity  scales 
are  based  on  damage  and  cause  effects,  it  follows  that  for  a  given  earthquake 


the  poorer  the  site  conditions  are,  the  larger  the  Intensity  will  be  at  a  given 
distance  from  the  source.  For  this  reason,  Cornell  and  Merz  (1975)  have  assigned  to 
the  parameter  C1  (In  Eq.  2.22)  values  of  2.6  for  rock  and  3.1  for  all  sites,  a  dif¬ 
ference  of  0.5  between  the  two  corresponding  intensities. 

Duration  of  earthquake  motion,  generally  speaking,  is  significantly  longer  on 
soil  sites  than  rock  sites,  as  observed  from  Table  2.10. 

B.  Near  Source  Effects 

Because  attenuation  laws  are  generated  from  data  observed  at  long  distances 
from  earthquake  sources,  the  validity  of  such  laws  regarding  short  distances 
from  earthquake  sources,  the  validity  of  such  laws  regarding  short  distances  from  a 
source  is  questionable.  Algermissen  and  Perkins  (1976)  suggest  that  acceleration 
tained  from  the  Schnabel  and  Seed  (1973)  curves. 

One  method  for  accounting  for  near  source  effects  on  the  law  is  the  use  of 
R  +  25  to  replace  R,  as  mentioned  earlier.  Donovan  and  Bornstein's  (1976)  recom¬ 
mendations  for  distance  dependent  parameters  b^ ,  and  b3  (Eq.  2.21)  are  another 
attempt  to  account  for  near  source  effects.  The  attenuation  law  suggested  by 
Milne  and  Davenport  (1969)  and  given  in  Table  2.5  implies  that  the  acceleration 
has  a  variation  of  e^'®^mR"2  at  long  distances,  and  of  e0’^11  at  short 
distances  from  the  source. 

Espinosa  (1977)  studied  the  intensity  attenuation  relationships  for  near 
and  intermediate  distances  using  the  data  from  the  1971  San  Fernando  earth¬ 
quake.  The  recomnended  relationship  Is  given  by: 

log  A  -  3.05  -  0.21  I  (2.30) 
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where  A  represents  the  epi central  distance  in  kilometres,  and  I  is  the  Modi¬ 
fied  Mercalli  Intensity.  Krinitzsky  and  Chang  (1977)  have  presented  relation¬ 
ships  between  acceleration  velocity  or  displacement  vs  MM  Intensity  for  near 
and  far  field  conditions. 

Cornell  and  Merz  (1975),  in  their  study  of  seismic  risk  in  Boston,  have 
assumed  that  there  is  no  attenuation  of  intensity  within  a  10-km  radius  around 
the  site.  Near  source  effects  on  acceleration  attenuation  are  not  yet  very 
well  understood.  Recent  papers  have  discussed  the  characteristics  of  near 
field  ground  motions  which  are  determined  based  on  physical  models  describing 
the  faulting  process  and  near  field  wave  propagation  (11,  15,  32,  52,  86,  99). 
Due  to  the  scarcity  of  near  field  data  for  large  magnitude  earthquakes,  reli¬ 
able  predictions  of  acceleration  close  to  the  source  cannot  yet  be  made. 

A  seismic  hazard  study  will  always  face  questions  regarding  the  proper 
selection  of  an  attenuation  law.  Following  careful  selection,  the  uncertain¬ 
ties  in  the  law  can  be  incorporated  in  the  hazard  analysis  following  the  pro¬ 
cedure  suggested  by  Cornell  and  Merz  (1975).  If  more  than  one  attenuation  law 
is  thought  to  be  applicable,  the  analysis  may  be  repeated  for  each  law  and  a 
Bayesian  estimate  of  seismic  hazard,  using  subjective  weights,  can  be  made. 

The  Bayesian  approach  will  be  described  in  the  following  section. 

Geologic  and  Geophysical  Information 

To  this  point,  discussion  has  been  limited  to  the  evaluation  of  local  and 
regional  seismicity,  based  on  parameters  obtained  from  statistical  data.  Prob¬ 
lems  and  limitations  involved  in  the  estimation  of  seismicity  based  on  seismo- 
logical  data  have  been  outlined.  In  sunwary,  existing  historic  records  are,  in 
most  cases,  too  recent  and  too  incomplete,  and  do  not,  therefore,  represent  the 
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true  nature  of  the  tectonic  process.  For  example,  the  absence  of  large  earth¬ 
quake  occurrence  indicated  in  the  historic  data  for  a  particular  region  may 
erroneously  be  Interpreted  as  an  indication  of  low  seismicity,  whereas  a  de¬ 
tailed  geologic  investigation  may  establish  the  existence  of  active  faults  in 
the  region  which  might  have  been  the  source  of  large  earthquakes  predating 
the  historic  record.  Thus,  studies  of  the  seismicity  of  a  region  based  soley 
on  statistical  Information  are  limited  in  the  conclusions  they  can  realistically 
or  safely  draw.  Bell  and  Hoffman  (1978)  have  made  a  comparative  study  between  the 
computed  seismic  hazard  based  on  geologic  evidence  and  the  computed  hazard  based 
on  Instrumental  data  for  a  site  near  Los  Angeles  (Fig.  2.20).  They  concluded  that, 
for  "normal  levels  of  acceptable  seismic  risk",  the  accelerations  which  were 
computed  based  on  geologic  evidence  exceeded  those  based  on  instrumental  data 
by  as  much  as  50%. 

It  is  now  accepted  that  geologic  and  geophysical  information  can  contribute 
significantly  to  the  determination  of  local  seismicity,  particularly  where  there 
exists  a  strong  correlation  between  past  earthquake  occurrence  and  the  surface 
faults  present,  such  as  in  California.  So  extensive  is  the  information  regard¬ 
ing  faults  and  fault  activity  that  Allen  (1976)  suggests  that  a  seismic  zoning 
map  for  California  based  on  geologic  data  alone  would  be  more  reliable  than  a 
map  based  on  purely  statistical  data.  However,  the  most  rational  procedure  in 
seismic  hazard  analysis  is  to  utilize  the  various  types  of  data  stemming  from 
different  sources.  Incorporating  the  subjective  judgments  of  the  geologist  and 
the  seismologist.  Baye's  theorem  provides  a  rational  procedure  toward  this  end. 

Seismic  hazard  analysis  following  the  Bayesian  approach  was  first  intro¬ 
duced  by  Esteva  (1969).  Various  other  researchers  have  since  discussed  the  ap¬ 
plication  of  Baye's  theorem  to  the  evaluation  of  local  seismicity.  However, 
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NUMBER  OF  OCCURRENCES  EXCEEDING  V 


only  recently  have  complete  formulations  of  the  Bayesian  approach  and  its  practi¬ 
cal  applications  In  seismic  hazard  analysis  been  properly  and  extensively  pre¬ 
sented. 

In  brief,  the  Bayesian  approach  requires  the  cooperative  effort  and  input 
from  geologists,  geophysicists  and  earthquake  engineers.  The  pertinent  geo¬ 
logic  Information  includes  the  Identification  of  all  active  faults  in  the  given 
region  and  the  assignment  of  activity  levels  for  each  fault.  Such  an  evaluation 
of  fault  activity  Involves  various  Investigative  approaches,  including  geologic 
mapping,  interpretation  of  aerial  photographs  and  remote  sensing  imagery,  direct 
subsurface  investigation,  geophysical  surveys,  age-dating  of  geologic  material, 
geodetic  measurements,  and  seismological  studies  (22,  59). 

Various  definitions  of  fault  activity  are  present  in  the  literature  (22,  59). 
Table  2.13  presents  the  fault  classification  used  by  Bell  and  Hoffman  (1978).  In 


TABLE  2.13 


FAULT  ACTIVITY  CLASSIFICATION  (After  Bell  and  Hoffman  1978) 


Fault  Activity 

Classification  Geologic  Time 


Historically  active 
Active 

Potentially  Active 
Inactive 


Historic 
Hoi ocene 
Pleistocene 
Pre-Pleistocene 


Displacement  Age 
(Years  before  Present! 


0  to  200 
200  to  11,000 
11,000  to  3,000,000 
Older  than  3,000,000 


addition  to  the  classification  of  a  fault,  the  degree  of  activity  for  each  fault 
must  be  determined,  most  often  In  terms  of  slip  rate,  cumulative  slip  during  a 
given  interval,  and  amount  of  slip  during  one  event  (22). 

The  recurrence  intervals  of  given  magnitude  events  can  be  determined  by 
the  estimation  of  creep  and  energy  liberated  by  shocks,  based  on  measured  fault 
displacements  per  unit  time  in  the  recent  geologic  time  scale  (43,  44).  Slemmons 
(1977)  has  presented  recurrence  Intervals  as  a  function  of  earthquake  magnitude 
and  strain  rates  across  fault  zones. 
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Geologic  data  for  a  region  can  also  provide  a  basis  for  estimating  the  max 


Imum  credible  earthquake,  m^ ,  that  can  occur  along  a  known  fault.  Such  an  esti¬ 
mate  is  usually  based  on  the  length  of  the  fault  rupture  and  the  type  of  fault 
movement.  Various  relationships  between  fault  rupture  length  and  earthquake 
magnitude  have  been  sumnarized  by  Slemmons  (1977).  Cluff  (1978)  warns  that  es¬ 
timates  of  m.j  made  soley  on  the  basis  of  the  length  of  fault  rupture  may  be  un¬ 
realistic,  and  recommends  consideration  of  the  following  five  factors  in  order 
to  obtain  more  reliable  estimates: 

1)  Geologic  evaluation  of  the  regional  tectonic  framework. 

2)  History  of  seismic  activity  along  the  fault  and  the  surrounding  region. 

3)  Geologic  history  of  displacement  along  the  fault. 

4)  Relationship  between  earthquake  magnitude  and  fault  rupture  length. 

5)  The  relationship  between  earthquake  magnitude  and  the  amount  of  fault 
displacement. 

Bayesian  Approach 

Due  to  variations  and  uncertainties  in  the  seismic  parameters  x(m),  B,  and 
m^,  estimated  on  the  basis  of  statistical  data,  these  parameters  are  often  treated 
as  random  variables  (16,  41,  43,  81,  82,  102).  Reliable  probability  distribu¬ 
tions  for  these  parameters  can  be  determined  by  incorporating  geologic,  geophys¬ 
ical  and  statistical  data,  along  with  professional  judgment,  into  a  cohesive  in¬ 
formation  source.  Application  of  Baye's  theorem  allows  for  such  a  rational  com¬ 
bination  of  the  widely  varying  data  sources  (10,  89). 

For  illustrative  purposes,  the  parameter  x(m),  defining  the  annual  rate  of 
occurrence  of  earthquakes  with  magnitudes  equal  to  or  greater  than  m  will  be 
used  to  present  the  basic  Bayesian  probability  theory  applied  in  seismic  hazard 
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analysis.  In  Eq.  2.2,  an  exponential  function  was  used  to  relate  x(m)  to  m. 

The  application  of  Baye's  theorem  requires  that  x(m)  be  treated  as  a  random 
variable  and  not  a  function  of  m.  This  can  be  accomplished  by  assuming  that 
the  x's  for  various  magnitudes  are  independent  of  one  another.  Thus,  Baye's 
theorem  can  be  systematically  applied  for  each  X(m^),  independently  for  various 
values  of  m^.  The  updated  values  of  x(mi)  obtained  by  the  application  of  Baye’s 
theorem  will  indicate  discontinuity  when  plotted  against  m^.  This  can  be  cor¬ 
rected  by  fitting  a  continuous  curve  through  the  computed  x(m^)  values  (82). 

For  purposes  of  convenience,  let  x(m.j)  be  denoted  by  x^ .  Baye's  theorem 
can  be  expressed  in  mathematical  terms  as  in  Eq.  2.31: 

f(data|x.)  f(x.) 

f(x  (data)  =  - - - ^ -  (2.31) 

J  f(data|xi)  f(x.)  dxi 

where  f(x^)  represents  the  initial  or  prior  estimate  of  the  probability  density 
function  of  the  random  variable  x^;  f(data|x^)  represents  the  sample  likelihood 
function  defining  the  probability  of  observing  the  data,  given  x^;  and  f ( x . | data ) 
is  the  posterior  or  improved  probability  density  function  of  x.,  given  the  data. 
The  numerator  on  the  right  side  of  the  equation  is  a  normalizing  factor,  so  that 
when  f(x^)data)  is  integrated  over  all  x^,  the  summation  will  equal  1.0. 

Eq.  2.31  is  applied  in  the  following  manner.  A  prior  distribution  for  x^ 
is  chosen.  The  choice  of  this  distribution  is  difficult  and  at  times  even  arbi¬ 
trary,  yet  Rosenblueth  (1969)  has  stated  in  a  discussion  of  the  Bayesian  approach 
that  this  difficulty  "is  a  poor  reason  to  drop  this  approach  in  favor  of  irra¬ 
tional,  even  if  traditional,  statistical  methods."  The  prior  distribution  most 


commonly  chosen  for  X,.  Is  ganrna  with  parameters  X^  and  K.j,  obtained  from  either 
geologic  or  subjective  input.  Newnark  and  Rosenblueth  (1971)  discuss  procedures 
for  evaluating  such  prior  information  based  on  geologic  and  seismological  data 
corresponding  to  the  three  macrozones  shown  in  Fig.  2.7. 

The  sample  likelihood  function  can  be  obtained  by  using  statistical  data 
on  earthquake  occurrence,  and  assuming  that  earthquake  occurrence  generally  fol¬ 
lows  a  poisson  process.  If  the  available  data  for  any  given  source  indicates 
that  during  a  period  of  time  t  years  in  length,  n  number  of  earthquakes 
of  magnitude  equal  to  or  greater  than  m  have  occurred,  then: 

(x.t)n  e"*** 

f(data|x.)  =  — 1 -  (2.32) 

1  nl 


Substituting  Eq.  2.32  and  the  prior  probability  with  the  parameters  x^  and 
K.j  in  Eq.  2.31,  and  evaluating  the  normalizing  factor,  the  posterior  probability 
density  function,  f(x..|data),  may  be  evaluated.  Mortgart  et  al.  (1977)  present 
results  of  the  calculations  which  indicate  that  the  posterior  density  function  of 
x..  is  also  gamma  with  parameters 
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where  Xg  and  Kg  are  the  posterior,  "updated"  values  of  the  parameters  x1  and 
Kj ;  and  y  and  are  the  mean  and  variance,  respectively,  of  X^. 

Newmark  and  Rosenblueth  (1971)  explain  the  application  of  Baye's  theorem 
where  the  prior  information  on  X..  is  given  as  n-j  expected  occurrences  of  mag¬ 
nitudes  of  m  >  m.  in  a  period  of  t^  years.  Assuming  a  gamma  distribution  for 

-1  /2 

X..  and  a  prior  coefficient  of  variation  of  X^  equal  to  n^  ,  the  expected 
value  and  the  variance  of  the  posterior  density  function  for  x^  are  given  by: 


u 


+  n 


+  t 


2  nl  +  n 

a  =  -  9 

(t,  +  tr 


(2.34) 


Making  the  assumption  that  the  prior  coefficient  of  variation  of  X^  is  equal 
-1/2 

to  n^  '  is,  in  effect,  equivalent  to  assigning  to  K-j  a  value  of  n^ .  Thus, 

Eqs.  2.33  and  2.34  yield  the  same  result.  Eq.  2.33  expresses  the  posterior 

results  in  general  terms  as  a  function  of  the  prior  parameter  for  the  gamma 

distribution,  whereas  Eq.  2.34  is  valid  only  if  the  prior  coefficient  of  the 

-1/2 

variation  of  x^  is  assumed  to  be  equal  to  n^  ,  or  essentially,  if  K^  is 
equal  to  n^. 

Thus,  application  of  Baye's  theorem,  which  utilizes  statistical  data,  im¬ 
proves  the  initial  or  prior  estimation  of  x.,  which  is  based  on  conceptual 
models  and  non-statistical  information.  As  resources  of  statistical  data  in¬ 
crease,  Baye's  theorem  can  be  applied  to  continuously  update  the  seismicity 
evaluation  of  a  given  region. 
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In  general,  the  value  of  the  contribution  of  statistical  information  to 
the  Bayesian  estimation  is  dependent  on  the  extent  of  the  available  historic 
data,  as  well  as  on  the  degree  of  uncertainty  in  the  prior  estimate.  If  the 
available  seismological  records  are  both  relatively  long  and  reliable,  the 
posterior  estimates  should  not  differ  too  greatly  from  those  based  on  the  sta¬ 
tistical  data.  However,  if  the  prior  estimates  are  based  on  thorough  assess¬ 
ments  of  geologic  evidence,  then  the  statistical  data  should  not  significantly 
change  the  prior  estimates  (42). 

In  Eq.  2.1  the  estimated  probability  of  n  events  with  magnitude  m  '  m^ 
occurring  in  t  years  was  conditional  on  X^,  and  because  X^  is  a  random  var¬ 
iable  with  a  distribution  of  f(\^|data),  the  unconditional  Bayesian  estimate 
of  this  probability  will  be  given  by  Eq.  2.35: 

P(n)  *  J*  P(n,  Xj)  f(Xi|data)dXi  (2.35) 

Bayesian  Estimate  of  Seismic  Hazard 

The  presence  of  uncertainties  in  several  of  the  parameters  used  and  in  the 
steps  followed  in  a  seismic  hazard  analysis  have  been  emphasized  in  this  Report. 
Cornell  and  Merz  (1975)  have  employed  a  rather  simple  method  to  account  for  these 
uncertainties  and  to  allow  for  Incorporation  of  subjective  evaluation  in  the 
hazard  studies.  The  procedure  involves  repeating  the  analysis  for  a  number  of 
alternative  values  of  the  questionable  parameters  and  assumptions.  The  result¬ 
ant  hazard  estimate  is  made  by  assigning  a  subjective  weight  to  each  alternative 
and  sunming  the  weighted  results  for  all  the  alternatives.  This  estimate  is 
referred  to  as  a  Bayesian  estimate,  which  Cornell  and  Merz  state  "reflects  the 


uncertainty  induced  both  by  the  behavior  of  nature  and  by  man's  incomplete 
state  of  knowledge  about  that  behavior."  Fig.  2.2 1  shows  the  subjective 
weights  assigned  by  Cornell  and  Merz  (1975)  to  the  various  assumptions  made 
in  their  study  of  seismic  risk  in  Boston.  In  this  example,  three  assumptions 
were  made  on  the  attenuation  law,  three  on  source  configuration,  and  four  on 
the  upper  bound  magnitude.  The  relative  weights  for  each  of  the  36  studies 
are  calculated  by  multiplying  the  subjective  weights  corresponding  to  each 
outcome.  For  example,  for  the  most  likely  hazard  curve,  the  relative  weight 
is  equal  to  0.5  x  0.5  x  0.3  =  0.075.  The  resulting  Bayesian  weighted  estimate 
curve  is  then  obtained  by  sunning  the  products  of  each  of  the  36  hazard  curves 
by  the  relative  weight  associated  with  that  curve. 

Source  Upper  Bound 

Attenuation  Geometry  Intensity  Relative 

Alternatives  Configurations  Perturbation  Weights 

0.075 
0.075 
0.05 
0.05 


FIG.  2.21  Subjective  Relative  Weights  Corresponding  to  Differ¬ 
ent  Alternatives  Used  by  Cornell  and  Merz  (1975) 

Fig.  2.22  shows  a  comparison  between  the  most  likely  estimate  of  risk  and 
the  Bayesian  estimate.  The  two  results  differ  significantly  at  site  intensities 
greater  than  VII,  a  discrepancy  which  is  attributed  to  the  perturbation  in  which 
the  upper  bound  intensities  in  the  near  source  areas  were  raised  from  (7.7,  8.7) 
to  XII. 
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ANNUAL  RISK  OF  EQUALING  OR  EXCEEDING  INTENSITY 


FIG.  2.22  Comparison  Between  Most  Likely  Estimate 
of  Risk  and  Bayesian  Estimate 
(From  Cornell  and  Merz  1975) 
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PART  3:  SEISMIC  HAZARD  CALCULATIONS 


The  calculated  seismic  hazard  at  a  site  can  be  expressed  in  terms  of  a 
curve  which  represents  the  probability  of  exceeding  a  given  seismic  parameter 
at  the  site,  vs  the  parameter  itself.  Previous  discussion  has  provided  a  pro¬ 
cedural  basis  for  the  performance  of  the  hazard  analysis.  Fig.  3.1,  which  is 
reproduced  after  Shah  (1978),  gives  a  schematic  outline  of  the  major  steps  in¬ 
volved  in  hazard  calculations.  In  Step  1,  potential  earthquake  sources  are 
identified,  and  the  length  or  area,  position  and  location  of  each  source  are 
established.  Step  2  involves  the  determination  of  the  recurrence  relation¬ 
ship  for  each  source,  based  on  geologic,  geophysical  and  historic  data.  In 
Step  3,  an  attenuation  law  is  selected  in  order  to  describe  the  seismic  para¬ 
meter  thought  to  best  represent  the  ground  motion  intensity  at  the  site.  Fin¬ 
ally,  Step  4  involves  the  actual  calculations  of  the  seismic  hazard  and  the 
estimation  of  the  probability  of  exceedence  (or  non-exceedence)  within  a  given 
time  period  of  various  values  of  the  selected  parameter. 

The  next  section  of  the  Report  will  present  a  brief  discussion  regarding  a 
number  of  the  basic  probability  concepts  used  in  seismic  hazard  studies,  fol¬ 
lowed  by  a  presentation  of  the  methodology  and  the  mathematical  formulations 
used  in  Step  4  of  the  analysis. 

Basic  Probability  Concepts 

The  recurrence  relationship  given  in  Step  2  of  Fig.  3.1  is  established  by 
a  study  of  the  local  seismicity.  This  relationship  was  expressed  in  Part  2 
in  terms  of  two  parameters:  1)  x(mQ),  the  annual  number  of  earthquake  occur¬ 
rences  with  magnitudes  equal  to  or  exceeding  mQ,  representing  the  minimum  mag- 
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nitude  of  engineering  interest,  and  2)  e,  the  magnitude-frequency  parameter 
which  defines  the  conditional  probability  of  exceedence,  P[M  >  m|E],  as  given 
by  Eq.  2.5.  Seismic  hazard  calculations  are  performed  in  order  to  obtain  an 
estimate  of  the  probability  that  at  least  one  event  will  occur  in  t  years, 
causing  the  value  of  a  parameter  A  to  exceed  a.  This  estimated  probability 
figure  is  arrived  at  by  calculating  the  conditional  probability  of  A  exceeding 
a,  given  that  an  earthquake  E  will  occur: 

P[A  >  a|E]  =  P[M  >  m|E]  (3.1) 

where  m  is  determined  from  the  attenuation  law  (which  relates  A  to  M  and  R), 
using  A  =  a.  The  conditional  probability  shown  in  the  right-hand  side  of  Eq. 

3.1  can  be  evaluated  using  Eq.  2.5,  and  thus  the  value  of  P[A  >  a|E]  is  estimated. 

If  *(m0)  earthquakes  are  expected  to  occur  annually,  the  expected  number  of 
special  events  causing  A  to  exceed  a  is  given  by: 

X(a)  =  x(mQ)  P[A  >  a|E]  (3.2) 

Provided  that  earthquake  occurrence  is  assumed  to  follow  a  poisson  process 
of  arrival,  these  special  events  with  an  annual  rate  of  occurrence  of  x (a )  will 
also  follow  a  poisson  process,  assuming  that  they  are  independent  events.  There¬ 
fore,  the  unconditional  probability  of  at  least  one  event  causing  A  >  a  occurring 
within  a  period  of  t  years  is  given  by: 

P[A  >  a,t]  at  least  one  =  1  -  e"x^a^  (3.3) 
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The  annual  probability  of  exceedence  will  then  be  given  by: 


P[A  >  a,  t  -  1]  =  1  -  e"X(a)  (3-4) 

-  A(a);  for  small  x(a) 

If  the  annual  probability  of  exceedence  is  known,  then  the  corresponding 
probability  for  t  years  can  be  obtained  following  one  of  two  alternate  pro¬ 
cedures.  One  method  is  to  calculate  A(a)  from  Eq.  3.4  and  simply  substitute 
it  in  Eq.  3.3  to  obtain  the  probability  for  t  years.  Alternatively,  P[A  >  a,t] 
can  be  calculated  directly  from  the  annual  probability,  using  Eq.  3.5: 

P[A  >  a,t]  =  1  -  (1  -P[A  >  a,t  =  l])1  (3.5) 

=  P[A  >  a,t  =  1]  x  t;  for  probabilities  less  than  0.1 

These  two  procedures  are  essentially  the  same  and  should  yield  identical  results. 

In  the  literature,  the  estimated  probability  of  a  seismic  event  occurring 
is  often  expressed  in  terms  of  its  return  period.  The  return  period,  Tr(a),  of 
an  event  causing  A  to  exceed  a,  is  defined  as  the  inverse  of  the  annual  probability  of 
at  least  one  event  occurring  (10): 


for  small  A (a) 


(3.6) 
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A  second  parameter  used  In  seismic  hazard  studies  involves  the  average  in¬ 
ter-arrival  times,  referred  to  as  the  recurrence  interval,  T  (a),  which  is  de- 

d 

fined  by  the  inverse  of  the  mean  number  of  annual  events  causing  A  to  exceed  a: 

T  (a)  =  - L_  (3.7) 

3  X(a) 

These  two  parameters  are  commonly  used  interchangeably,  and  therefore 
have  been  confused.  A  comparison  between  Eq.  3.6  and  Eq.  3.7  leads  to  the 
conclusion  that  the  return  period  and  the  recurrence  interval  are  identical 
only  if  the  computed  probability  is  low.  Fig.  3.2  gives  the  ratio  of  return 
period  to  recurrence  interval  as  a  function  of  the  mean  annual  number  of  events, 
x .  It  is  clear  that  the  return  period  and  the  recurrence  interval  are  compar¬ 
able  only  if  x  is  less  than  0.1,  or  T  (a)  is  greater  than  10  years.  In  addi- 

a 

tion,  the  concept  of  recurrence  interval  is  often  misinterpreted  (2).  For  ex¬ 
ample,  the  probability  of  at  least  one  event  with  a  recurrence  interval  of 
T  (a)  occurring  in  T  (a)  years  is  not  equal  to  1.0,  as  might  be  expected,  but 

O  O 

rather  is  equal  to  0.63,  as  shown  in  Eq.  3.8: 

P[A  >  a,  t  =  T  (a)]  =  1  -  e'X(a)  *  1/x(a)  =  0.63  (3.8) 

a 

In  order  to  avoid  confusion,  it  is  advised  that  the  computed  hazard  results  be 
stated  in  terms  of  the  probability  of  exceedence  within  a  given  period  of  time. 
This  manner  of  presenting  the  results  provides  for  a  direct  association  of  the 
hazard  to  the  design  life  span  of  a  structure. 
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RECURRENCE  INTERVAL,  Ta,  YEARS 


Mathematical  Formulations 


A  seismic  hazard  analysis  can  employ  three  general  earthquake  source  types 
in  order  to  describe  the  potential  future  seismic  activity  for  a  region.  This 
section  of  the  Report  presents  the  basic  methodology  involved  in  the  analysis 
considering  the  three  source  types,  which  are  referred  to  as  point,  line  and 
area  sources.  For  reasons  of  simplicity,  the  attenuation  law  will  be  assumed 
to  be  deterministic,  having  the  form  A  =  b^e^2m(R  +  25)”b3,  and  it  will  be  further 
assumed  that  m^  =  ».  A  general  procedure  for  a  hazard  analysis  which  employs  a 
probabilistic  rather  than  deterministic  attenuation  function  is  presented  by 
McGuire  (1976)  and  summarized  in  Appendix  A. 

Point  Source 

Fig.  3.3  shows  a  site  located  R  km  from  a  point  source.  The  conditional 
probability,  P[A  >  a|E]  of  acceleration  A  exceeding  an  assumed  value,  a,  given 
that  an  earthquake  E  will  occur,  is  obtained  by  Eq.  3.9: 

P[A  >  a|E]  =  P[M  >  m|E]  (3.9) 

where  a  >  b^eb2mO(R  +  25)~b3 
*  -h 

and  m  =  ln[  t-  (R  +  25)  °3] 


Employing  the  complementary  cumulative  distribution  function  for  m,  given 
by  Eq.  2.5,  the  conditional  probability  of  A  exceeding  a  can  be  expressed  by: 


P[A  >  a|E] 


exp  J 


K  -1 /bp 

•B[ln[4-  (R  +  25)'b3] 


(3.10) 
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IINE  SOURCE 


SITE 

TOP  VIEW 


FIG.  3.3  Definitions  of  Source  Geometry  Parameters 
(After  Shah  et  al.  1975) 


=  eBmo 


(  )-e/b2  (R  +  25)‘b3e/t>2 


Therefore,  the  mean  annual  number  of  seismic  events  causing  acceleration 
A  to  exceed  a  is  given  by  Eq.  3.11: 

A(a)  =  x(m0)  [  eBmo  (  -i-  )_B/b2  (R  +  25)‘b3e/b2  ]  (3. 11) 

where  x(m0)  represents  the  mean  annual  number  of  earthquakes  with  M  >  mQ.  As¬ 
suming  a  poisson  process  of  arrival  for  these  events,  the  annual  probability 
of  occurrence  at  the  site  of  at  least  one  event  having  A  >  a  is  given  by: 

P[A  >  a,  t  =  1]  =  1  -  exp[-A(mJeBm°(  -4-  )'B/b2  (R  +  25)"b36/b2]  (3.12) 

0  D-j 

Line  Source 

The  line  source  was  originally  proposed  as  a  model  for  a  potential  earth¬ 
quake  source  by  Cornell  (1968)  and  is  employed  in  cases  where  the  epicenters  of 
historic  earthquakes  fall  along  a  line,  or  when  geologic  and  geophysical  inves¬ 
tigations  indicate  the  presence  of  an  active  fault.  Shah  (1975)  describes  one 
method  for  modeling  known  faults,  which  divides  the  line  representing  the  fault 
into  k  small  segments  of  length  dL,  as  shown  in  Fig.  3.3.  Each  segment  is 
then  treated  as  a  point  source,  with  a  hypocentral  distance  equal  to  the  dis¬ 
tance  of  the  segment  from  the  site.  An  attenuation  law  is  then  applied,  re¬ 
lating  acceleration  to  earthquake  magnitude  and  hypocentral  distance,  and  the 
calculations  for  the  hazard  analysis  are  as  follows.  Based  on  the  assumption 
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that  a  future  earthquake,  should  it  occur,  has  an  equal  likelihood  of  occur¬ 
rence  anywhere  along  the  entire  length  of  the  fault,  the  mean  annual  rcte  of 


earthquake  occurrence  for  each  segment,  i,  will  be  equal  to  x(mo)dl,  where 
x(mo)  represents  the  mean  annual  rate  of  occurrence  of  earthquakes  having 
M  >  mQ  per  unit  length. 

The  contribution  of  each  segment  i  to  the  mean  number  of  events  with 
A  exceeding  a  can  be  expressed  by: 


^(a)  =  X(mo)  dl  [eBm°  (  )’B/b2  (r  +  25)'b36/b2  ] 


(3.13) 


From  Fig.  3.3,  R  can  be  expressed  in  terms  of  the  perpendicular  distance 
to  the  fault,  the  focal  depth  and  L^,  the  distance  of  the  segment  i  from  the 
perpendicular  on  the  line  source: 


R  ■  (02  *  L(2  *  H2) 


(3.14) 


Substituting  the  value  of  R  into  Eq.  3.13  and  adding  the  contributions  of 
all  the  segments  to  the  mean  rate  of  occurrence,  x(a), 

X(a)  =  X(m0)  eemo  (  -jj-  )‘6/b2  E  [(D2  +  L*  +  H2)1/2  +  25]"b3e/b2  dLi  (3 


Thus,  the  annual  probability  of  A  exceeding  a  is  given  by: 


P[A  >  a,  t  =  1]  =  1  - 


-x(a) 


(3.16) 
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Arguments  have  been  presented  that  such  a  treatment  of  a  fault,  i.e.,  divid¬ 
ing  it  into  small  segments  or  into  areas  as  described  in  Part  2,  does  not  con¬ 
sider  the  fault  as  a  line  source.  Rather,  dividing  a  fault  into  small  segments 
and  using  an  attenuation  law  which  is  based  on  hypocentral  distance  (in  this 
case,  distance  of  the  site  from  the  segment)  is  actually  equivalent  to  treat¬ 
ing  the  fault  as  a  series  of  point  sources.  For  this  reason,  as  well  as  for 
computational  efficiency,  Donovan  (1978)  has  preferred  to  use  only  point  sources 
in  seismic  hazard  analysis. 

DerKiureghian  and  Ang  (1977)  have  described  the  problems  associated  with 
modeling  a  fault  as  a  point  source.  Such  so-called  "point  source"  models 
for  faults  assume  that  during  an  earthquake,  energy  is  radiated  from  the  focus 
of  the  earthquake.  While  such  an  assumption  is  valid  for  small  earthquakes, 
during  a  large  magnitude  earthquake  the  energy  released  would  be  distributed 
along  the  rupture  zone  which  may  be  propagated  close  to  the  site  under  study. 
Indicating  that  ignoring  the  effect  of  the  rupture  length  would  tend  to  under¬ 
estimate  the  real  hazard,  Der  Kiuregfian  and  Ang  (1977)  have  proposed  a  re¬ 
search  model  for  faults  which  is  based  on  the  assumption  "That  an  earthquake 
originates  at  the  focus  and  propagates  as  an  intermittent  series  of  fault  rup¬ 
tures  or  slips,  and  that  the  maximum  intensity  of  ground  shaking  at  a  site  is 
determined  by  the  slip  that  is  closest  to  the  site." 

DerKiureghian  and  Ang  (1977)  have  outlined  three  types  of  sources  in  de¬ 
scribing  the  fault-rupture  model,  which  are:  1)  a  well-defined  fault  line; 

2)  an  area  of  known  fault  direction;  and  3)  an  area  of  unknown  faults.  A 
brief  presentation  of  the  type  1  fault-rupture  model  follows. 

Fig.  3.4  shows  the  location,  position  and  length  of  a  known  fault.  The 
distance  x  of  the  focus  from  the  perpendicular  line  to  the  fault  line  is 
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FIG.  3.4  Definitions  of  Parameters  for  the  Fault-Rupture  Model 
(From  DerKlureghlan  and  Ang  1977) 


given  by  Eq.  3.17: 


x  -  (r2  -  r02)  1/2  +  s/2  (3.17) 

where  s  represents  the  length  of  the  slip,  which  Is  related  to  earthquake  mag¬ 
nitude: 


s  *  exp(am  -  b) 


(3.18) 


where  a  and  b  are  constants  (97).  If  It  Is  assumed  that  the  focus  of  a  future 
earthquake  would  have  an  equal  likelihood  of  occurrence  anywhere  along  the  length 
of  the  fault,  the  probability  distribution  fR(r)  for  r  ,  the  minimum  distance 
from  the  site  to  the  slip,  can  be  derived  from  Eq.  3.17  and  the  probability  dis¬ 
tribution  for  x 


fx(x) 


j  0  S  X  S  l 


otherwise 


(3.19) 


After  the  distribution  fR( r)  has  been  determined,  the  conditional  proba¬ 
bility  distribution  for  A  ,  site  Intensity,  given  that  an  earthquake  of  mag¬ 
nitude  m  occurs  In  source  1  ,  can  be  obtained  from  the  attenuation  function 
which  relates  A  to  m  and  r  .  The  probability  that  A  exceeds  a  ,  con¬ 
ditional  on  the  occurrence  of  an  earthquake  at  source  1  ,  Is  obtained  by  in¬ 
tegration  over  the  appropriate  values  of  a  and  m  : 


P  [A  >  a  |  E] 


,M1  max 

/  / 


f(a)  fM(m)  dadm 
A  |  E,m  M 


(3.20) 


The  annual  probability  of  exceedence  can  then  be  obtained  using  Eqs. 


3.2  and  3.4. 


Area  Source 


If  In  a  particular  region  the  epicenters  of  historic  earthquakes  are  scat¬ 
tered  randomly  due  to  the  presence  of  numerous  faults  and/or  errors  In  the  es¬ 
timations  of  eplcentral  locations  of  past  earthquakes,  an  area  source  model  may 
be  most  appropriate  to  describe  the  future  seismicity  of  the  region.  If  the 
rate  of  occurrence  per  unit  area,  the  focal  depth  and  the  attenuation  are  the 
same  for  all  earthquakes  originating  within  an  area  source,  that  source  is  com¬ 
monly  referred  to  as  a  "uniform"  source.  For  computational  efficiency,  an  area 
source  is  often  represented  by  a  segment  of  an  annular  ring,  as  shown  in  Fig. 

3.5.  If  the  area  source  considered  In  an  analysis  does  not  have  an  annular 
shape,  the  source  is  divided  into  subsources  having  annular  shapes,  as  shown 
in  Fig.  3.6. 

The  annual  probability  calculations  for  an  area  source  are  as  follows. 

An  example  area  source  has  an  elemental  area  of  rdre,  with  a  hypocentral  dis- 
tance  of  (r  +  H£)'/2  (see  Fig.  3.5).  The  mean  annual  number  of  earthquakes 
occurring  in  this  area  is  given  by  A(mQ)  •  rdre,  and  the  conditional  probabil¬ 
ity  of  an  earthquake  occurring  with  A  exceeding  a  is  given  by  Eq.  3.10.  Therefore, 
the  annual  mean  number  of  events  with  A  exceeding  a  occurring  in  the  elemental  area 
Is  given  by: 

x(a)  -  X(mQ)  [eBm°  (  )'e/b2  (R  +  25)~b36/b2  ]  rdre  (3.21) 

The  total  mean  number  of  events  with  A  exceeding  a  occurring  in  the  annular  ring 
defined  by  the  two  radii  r-j  and  r ^  and  the  angle  e  is  given  by: 

x(a)  -  X(m  )eeBm<>  (  -4-  )'a/b2  \Z  [(r2  +  H2)1/2  +  25)"b36/b2  .  rdr  (3.22) 

1  rl 
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SOU*Cf 


TOP  VI  [W 


FIG.  3.6  Arc  Areas  Generated  by  Program  to  Assign  Seismicity 
Within  Quadrilateral  Source  Area  as  a  Function  of 
Oi stance  from  Site. 

(From  McGuire  1976) 


r*i  «-i 


The  estimated  annual  probability  of  A  exceeding  a  is  then  obtained  from: 


P[A  >  a]  =  1  -  e'X(a)  (3.23) 

The  total  annual  probability  at  a  site  is  then  obtained  by  summing  the 
contributions  of  all  sources  to  this  probability: 

N 

P[A  >  a]  =  1  -  e"i^1Xi(a)  (3.24) 

where  x^(a)  represents  the  annual  mean  number  of  events  with  A  exceeding  a  occurring 
due  to  source  i;  and  N  is  the  total  number  of  sources. 

When  the  uncertainty  in  the  attenuation  law  is  incorporated  in  the  analysis, 
the  integral  involved  in  the  calculation  of  the  annual  probabilities  can  be  cal¬ 
culated  numerically. 

The  incorporation  of  uncertainties  in  the  attenuation  law  results  in  more 
complex  mathematical  formulations.  In  Appendix  A,  McGuire's  (1976)  formulation 
of  the  mathematical  expression  used  in  such  a  case  is  briefly  presented. 

This  type  of  calculation  becomes  problematic  when  in  a  particular  study 
there  exists  a  large  number  of  sources  of  various  shapes  to  be  considered,  and 
the  uncertainites  of  the  various  parameters  are  to  be  incorporated  in  the  cal¬ 
culation.  Computer  programs  have  been  coded  for  the  purpose  of  performing 
these  calculations.  The  programs  most  widely  used  and  which  are  currently  avail¬ 
able  are:25 

1)  M.I.T.  program  by  Professor  C.  A.  Cornell,  Department  of  Civil  Engineer¬ 
ing,  Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139 

2)  USGS  program  by  Dr.  R.  K.  McGuire,  Branch  of  Earthquake  Tectonics, 

Denver  Federal  Center,  Denver,  Colorado  80225 
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PART  4:  APPLICATIONS 


This  Part  of  the  Report  will  discuss  the  various  applications  of  seis¬ 
mic  hazard  analysis.  A  review  of  the  literature  has  indicated  that  the  re¬ 
sults  of  seismic  hazard  studies  have  been  utilized  for  a  variety  of  purposes. 

Cornell  and  Merz  (1975)  conducted  a  seismic  hazard  analysis  in  Boston, 
the  result  of  which  became  a  major  contributing  factor  to  the  formulation  of 
a  new  Massachusetts  building  code  which  provides  for  earthquake  design.  Seis¬ 
mic  hazard  analysis  has  been  applied  by  many  investigators  throughout  the  world 
toward  the  development  of  hazard  or  zoning  maps  for  different  countries.  For 
the  United  States,  the  most  recent  hazard  maps  (presented  and  discussed  later) 
are  the  map  of  Algermissen  and  Perkins  (1976),  which  provides  probabilistic  es¬ 
timates  of  maximum  rock  acceleration,  and  the  map  developed  by  the  Applied  Tech¬ 
nology  Council  (1978)  providing  estimates  for  "effective  peak  acceleration" 
and  "effective  peak  velocity"  for  the  contiguous  48  states. 

Whitman  et  al.  (1975)  have  described  a  seismic  design  decision  analysis 
which  presents  a  methodology  for  evaluating  alternate  structural  design  re¬ 
quirements.  The  procedure  utilizes  the  results  from  a  hazard  analysis  expressed 
in  terms  of  Modified  Mercalli  Intensity. 

Current  earthquake  mitigation  efforts  and  earthquake  insurance  policy  for¬ 
mulations  consider  seismic  hazard  results  as  invaluable  input. 

Seismic  hazard  analysis  has  been  used  for  such  purposes  as  the  siting  of 
nuclear  power  plants  and  in  selecting  the  operating  basis  earthquake  (OBE)  (21, 
94).  Applications  to  other  types  of  facilities  Including  buildings,  dams  and 
offshore  facilities  are  also  reported  in  the  literature  (5,  6,  7,  8,  9,  17,  28, 
35,  65,  85,  104). 
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The  type  of  seismic  hazard  analysis  reviewed  herein  has  been  extended  to 
the  joint  probability  that  various  levels  of  earthquake  shaking  will  be  exper¬ 
ienced  at  a  number  of  discrete  sites  (98,  100,  106).  Such  modifications  allow 
for  the  study  of  seismic  risk  for  whole  communities  and  cities.  Special  atten¬ 
tion  has  been  focused  in  recent  years  on  the  seismic  risk  studies  of  lifelines 
(95).  Seismic  hazard  maps  are  most  valuable  in  the  risk  studies  of  networks 
of  various  lifelines  extending  over  large  seismically  active  regions. 

In  general,  a  seismic  hazard  analysis  provides  a  probabilistic  description 
of  the  seismic  load  on  a  facility  which  then  can  be  used  in  a  seismic  risk  analy¬ 
sis  using  a  probabilistic  resistance  function  to  arrive  at  a  probabilistic  estim¬ 
ate  of  damage  or  consequence  to  that  facility.  The  r’sk  analysis  can  be  described 
by  Eq.  4.1 : 


P[F]  -  J  P[F|P]  •  fp(p)dp  (4.1) 

where  fp(p)  is  the  probability  density  function  of  the  parameter  p.  This  para¬ 
meter  describes  the  hazard  at  the  site  and  constitutes  the  output  from  a  seis¬ 
mic  hazard  analysis.  P[F|p]  is  the  conditional  probability  of  failure,  damage 
or  consequence  on  the  uncertain  resistar-e  of  the  facility  to  the  hazard  p. 

This  probability  is  dependent  urt  ths  uncertain  resistance  of  the  facility  to  the 
hazard  p.  ?  total  urn.  .nditional  probability  of  failure  is  then  determined  by 

the  integral  of  the  product,  shown  in  Eq.  4.1,  over  all  values  of  interest  of 
the  hazard,  p. 

Very  often,  a  design  level  load  associated  with  an  acceptable  probability  of 
exceedence  is  chosen  from  the  output  of  a  seismic  hazard  analysis,  and  the  de¬ 
sign  of  a  facility  is  carried  out  in  a  conventional,  deterministic  way,  thus 
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avoiding  a  complete  risk  analysis  which  evaluates  the  contribution  to  this  risk 
of  a  range  of  values  of  the  input  load.  Donovan  (1978)  describes  the  common 
practice  of  selecting  two  levels  of  the  design  load  from  the  hazard  curves  for 
structural  design  purposes.  The  specific  seismic  event  for  which  a  structure 
is  to  be  designed  to  remain  elastic  should  that  event  occur  would  correspond 
to  a  50%  probability  of  exceedence  during  a  50-year  period.  The  adequacy  of 
the  design  in  terms  of  the  required  ductility  would  then  be  evaluated  using  the 
event  with  a  10%  probability  of  being  exceeded  in  a  50-year  period.  This  latter 
level  is  also  what  the  recent  hazard  maps  have  used. 

The  following  discussions  will  provide  descriptions  of  the  various  proce¬ 
dures  and  parameters  used  in  describing  the  output  of  the  seismic  hazard  analy¬ 
sis,  which  then  becomes  the  input  load  in  the  consideration  of  structural  and 
geotechnical  problems.  The  hazard  maps  of  Algermissen  and  Perkins  (1976)  and 
the  Applied  Technology  Council  (1978)  appear  separately  at  the  end  of  this  sec¬ 
tion. 


Input  to  Structural  Analysis 

This  Report  has  outlined  and  described  the  various  parameters  presently 
used  in  seismic  hazard  studies  to  describe  the  intensity  of  ground  motion  at  a 
site  during  a  seismic  event,  including  peak  ground  acceleration,  velocity,  dis¬ 
placement,  spectral  velocity,  Fourier  spectrum  ordinate,  and  duration.  Design 
considerations  for  a  particular  structure  may  require  input  from  one  or  more  of 
these  parameters.  In  most  cases,  peak  acceleration  and/or  velocity  is  not  suf¬ 
ficient  to  completely  describe  the  ground  motion  characteristic  (76).  Frequency 
content  and  duration  of  motion  are  significant  to  the  determination  of  the  dy¬ 
namic  response  of  a  structure.  Thus,  the  most  appropriate  and  desirable  seismic 


101 


parameter  is  the  spectral  velocity  for  the  structure,  or  preferably,  the  spec¬ 
tral  velocities  for  the  frequency  range  pertinent  to  the  structure.  There  may 
be  more  than  one  approach  followed  in  the  hazard  analysis  in  order  to  arrive 
at  this  design  spectrum,  depending  on  the  nature  of  the  available  seismological 
data  and  the  attenuation  law  chosen  for  the  analysis. 

The  available  seismological  data  is  generally  stated  in  terms  of  either 
Modified  Mercalli  Intensity  or  Richter  Magnitude.  Therefore,  there  will  normally 
be  two  approaches  to  the  estimation  of  the  probability  of  exceedence  of  Sv,  each 
of  which  will  be  discussed  separately. 

Intensity  Data 

Seismicity  data  which  is  stated  in  terms  of  epicentral  intensities  allows 
for  two  analytic  options. 

^  I0"n 

This  approach  necessitates  the  conversion  of  all  epicentral  data  to  earth¬ 
quake  magnitude  figures  using  the  relationships  given  in  Fig.  2.4,  and  the  analy¬ 
sis  is  then  proceeded  with  on  the  basis  of  earthquake  magnitudes.  A  more  de¬ 
tailed  discussion  of  this  procedure  will  be  presented  subsequently. 

2)  i  -  I  p  -  s 
1  *o  s _ _v 

This  approach  makes  use  of  epicentral  intensity  data.  The  seismic  hazard 
analysis  is  completed  expressing  the  results  in  terms  of  site  intensity.  The 
site  intensities  are  then  related  to  peak  acceleration,  velocity  or  displace¬ 
ment  by  means  of  the  published  relationships.  Krinitzsky  and  Chang  (1977)  have 
presented  such  relationships  for  near  and  far  field  conditions. 

Alternatively,  the  attenuation  law  selected  to  relate  epicentral  intensity 
to  site  intensity  can  be  combined  with  any  one  of  the  attenuation  laws  proposed 


by  McGuire  (1976)  (see  Table  2.5)  relating  site  intensity  to  acceleration.  A 
combined  relationship  can  then  be  developed  which  relates  epicentral  intensity 
directly  to  peak  ground  acceleration,  velocity  or  displacement.  The  output  from 
the  seismic  hazard  analysis  will  directly  yield  the  probability  of  exceedence  of 
a,  v  or  d.  The  question  of  which  parameter  (a,  v  or  d)  should  be  selected  for 
study  was  discussed  in  Part  2  of  the  Report.  In  brief  summary,  peak  ground  ac¬ 
celeration  may  be  used  if  the  period  of  the  structure  is  short,  while  peak  ground 
velocity  is  more  appropriate  for  intermediate  period  ranges  (1  cps).  Peak  ground 
displacement  is  the  recommended  parameter  for  use  in  design  considerations  for 
long  period  structures,  such  as  offshore  towers  (7,  29). 

In  cases  where  the  hazard  curve  is  expressed  in  terms  of  site  intensity, 
duration  of  ground  shaking  can  be  estimated  for  each  site  intensity  using  the 
relationship  given  in  Table  2.10.  The  estimation  of  the  design  response  spec¬ 
trum  is  made  by  scaling  the  standardized  response  spectra  using  the  design  level 
acceleration  and  velocity  corresponding  to  an  "acceptable"  probability  of  ex¬ 
ceedence  (or  non-exceedence)  within  a  given  time  period.  McGuire  (1977b)  has  in¬ 
dicated  that  such  spectra  "are  found  to  be  inconsistent  in  terms  of  risk  for 
building  sites  very  close  and  very  far  from  faults",  and  has  recommended  an 
alternate  procedure  which  will  be  discussed  subsequently. 

Magnitude  Data 

Two  alternative  approaches  to  hazard  analysis  are  available  when  the  seis¬ 
mic  data  for  a  region  is  stated  either  in  terms  of  earthquake  magnitude  or  in 
terms  of  epicentral  intensity,  converted  to  magnitude  for  the  purposes  of  the 
analysis. 
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This  approach  utilizes  a  magnitude-distance  attenuation  law  to  describe  the 
hazard  at  a  site  in  terms  of  a  parameter,  p,  which  can  be  either  peak  ground  ac¬ 
celeration,  velocity  or  displacement,  whichever  of  these  parameters  is  thought  to 
be  most  appropriate.  Considerations  pertinent  to  the  selection  and  use  of  these 
parameters  in  determining  a  design  response  spectrum  and  in  incorporating  site 
effects  were  previously  discussed.  McGuire  (1977a)  has  provided  a  function  re¬ 
lating  Fourier  amplitude  spectrum  to  earthquake  magnitude,  distance  and  soil  type 
presented  in  Part  2  in  Eq.  2.24.  A  Fourier  spectrum  gives  a  qualitative  measure 
of  the  nature  of  earthquake  motions.  The  Fourier  amplitude  spectrum  is  a  lower 
bound  to  the  pseudo-velocity  spectrum  for  zero  damping.  Employing  this  relation¬ 
ship  for  different  frequencies,  a  Fourier  amplitude  spectrum  associated  with  a 
chosen  fractile  can  be  developed. 

McGuire  (1977b)  has  also  proposed  a  parameter,  s,  which  is  proportional  to  Sv, 
for  use  in  scaling  normalized  spectra: 

s(l  Hz,  7%  damping)  =  0.155  x  10°*382m(R  +  25 )'0, 587 (in/sec)  (4.2) 

He  proposes  the  use  of  s  for  the  middle  frequency  range  (1  Hz)  and  the 
use  of  acceleration,  a,  for  the  high  frequency  range  to  scale  standardized 
Newmark-Blume-Kapur  response  spectra  (83).  Thus,  the  seismic  hazard  analysis 
is  conducted  using  attenuation  laws  for  s  and  a.  The  design  values  of  s  and 
a,  corresponding  to  an  "acceptable"  probability  level,  are  then  selected  from  the 
hazard  curve  and  applied  as  described  above.  In  addition,  design  values  of  s 
and  a  can  be  used  with  the  corresponding  attenuation  laws  chosen  for  the  analy- 


sis  to  define  a  design  earthquake  of  magnitude  and  hypocentral  distance  Rd< 
The  determination  of  a  design  earthquake  allows  for  the  estimation  of  the  dur¬ 
ation  of  ground  shaking  associated  with  the  design  criterion,  using  Figs.  2,21 
and  2,22. 


2)  m,  R  -*•  Sv 

A  seismic  hazard  analysis  following  this  approach  utilizes  an  attenuation 
law  in  order  to  relate  earthquake  magnitude  and  distance  to  Sy.  Thus,  a  direct 
estimation  of  the  probability  of  exceedence  of  Sy  for  a  particular  structure 
can  be  made  using  the  attenuation  constants  for  Sy  corresponding  to  the  period 
of  the  structure.  If  a  complete  response  spectrum  is  desired,  the  hazard  analy¬ 
sis  may  be  repeated  for  each  period,  changing  the  attenuation  constants  as  re¬ 
quired. 


While  the  parameter  Sy(T)  is  most  useful  in  structural  analyses,  most  soil 
dynamics  studies  in  earthquake  engineering  utilize  peak  ground  acceleration, 
duration,  or  a  time  history  of  ground  acceleration. 

It  is  well-known  that  the  susceptibility  of  soils  to  liquefaction  is  de- 
penderit  not  only  on  the  peak  ground  acceleration  (or  shear  stress),  but  also  on 
the  number  of  applications  of  the  acceleration  (or  stress).  Youd  and  Perkins 
(1978)  have  combined  cyclic  shear  strength  information  with  attenuation  laws 
for  cyclic  shear  stress  and  the  number  of  significant  cycles,  expressed  in  terms 
of  earthquake  magnitude  and  distance,  and  have  developed  a  criterion  for  lique¬ 
faction,  given  by  Eq.  4.3: 


m  =  C.j  +  C2  log  R 


(4.3) 
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where  C1  and  are  coefficients  which  vary  with  the  attenuation  factors  and 
local  site  conditions.  Eq.  4.3  can  be  applied  In  a  seismic  hazard  analysis  for 
a  site  such  that  for  each  segment  of  a  source  zone  having  a  distance  R^,  the 
annual  number  of  earthquakes  with  magnitudes  greater  than  Is  calculated. 

The  total  annual  mean  number  of  these  special  events  causing  liquefaction  at 
the  $1te  Is  obtained  by  sunmlng  the  contributions  to  this  number  of  all  seg¬ 
ments  from  each  source.  The  probability  of  n  such  occurrences  In  a  period  of 
t  years  can  then  be  obtained,  assuming  that  the  events  causing  liquefaction 
follow  a  polsson  process  of  arrival. 

McGuire  et  al.  (1978)  have  described  a  liquefaction  risk  analysis  In 
which  an  attenuation  law  for  the  earthquake- Induced  shear  stress  was  compared 
to  the  cyclic  shear  strength  of  the  soil  deposit  to  compute  the  probability 
of  liquefaction.  The  constants  for  the  shear  stress  attenuation  law  were  ob¬ 
tained  by  performing  a  regression  analysis  on  the  results  obtained  from  one¬ 
dimensional  wave  propagation  studies. 

Yegian  and  Whitman  (1978)  have  presented  a  risk  analysis  for  liquefaction 
which  employs  a  parameter  referred  to  as  the  liquefaction  potential  index,  LPI, 
given  by  Eq.  4.4: 

LPI  *  b1eb2m(R  +  25)~b3  (4.4) 

The  parameters  b-j,  b2  and  bj  are  obtained  based  on  the  soil  properties  at 
the  site  and  from  Interpertation  of  field  observations  of  cases  of  liquefaction 
and  no  liquefaction,  as  well  as  on  an  assumed  acceleration  attenuation 
law.  It  is  noted  that  the  form  of  LPI  Is  similar  to  the  form  of  the  attenuation 
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laws  commonly  used  in  seismic  hazard  analysis.  Since  the  LPI  is  lognormally 
distributed,  Eq.  4.4  can  be  used  as  an  attenuation  law,  and  the  probability  of 
liquefaction  (LPI  >  1.0)  can  be  calculated  employing  any  of  the  conventional 
procedures  of  seismic  hazard  analysis  reviewed  in  this  Report,  It  is  believed 
that  the  use  of  LPI  implicitly  takes  into  account  duration  of  earthquake  motion. 

When  the  duration  associated  with  a  design  level  acceleration  is  desired, 
the  seismic  hazard  analysis  is  often  repeated  using  an  attenuation  law  for  dur¬ 
ation.  In  reference  to  the  duration  considerations  discussed  previously,  it 
must  be  brought  to  attention  here  that  separate  consideration  of  duration  esti¬ 
mates  falsely  assumes  that  acceleration  and  duration  are  independent  variables. 
Duration  and  peak  ground  acceleration  are  both  related  to  earthquake  magnitude 
and  distance,  the  parameters  used  to  describe  the  essential  characteristic  of 
an  earthquake;  and,  therefore,  they  are  not  statistically  independent  variables. 
The  procedure  proposed  by  McGuire  (1977b)  and  described  earlier  regarding  the  cal 
culation  of  a  design  earthquake  with  magnitude  md  and  distance  Rd  can  be  used 
to  relate  duration  to  design  acceleration  more  appropriately,  using  the 
relationships  presented  by  Chang  and  Krinitzsky  (1977). 

In  general,  soil  dynamics  problems  associated  with  major  structures,  includ 
ing  power  plants  and  dams,  require  the  use  of  a  time-history  record  of  motion 
on  firm  ground.  The  proper  selection  of  a  time-history  consistent  with  a  spe¬ 
cific  design  criterion  is  difficult.  There  presently  exists  no  systematic  pro¬ 
cedure  for  this  selection;  in  lieu  of  a  tested  procedure,  all  pertinent  infor¬ 
mation,  including  available  instrumental  records  of  major  aftershocks,  should  be 
weighed  with  careful  professional  judgment  in  order  to  select  a  time-history. 

In  most  cases,  more  than  one  time-history  is  selected  and  the  influence  of  the 
variation  in  the  time-history  characteristics  is  evaluated.  Among  the  para- 
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meters  used  in  the  selection  process  of  a  time-history  are  peak  acceleration  and 
the  ratio  of  peak  acceleration  to  peak  velocity.  McGuire  (1978)  has  pointed 
out  that  this  ratio  is  not  a  constant,  but  rather  is  dependent  on  the  distance 
of  the  site  from  a  potential  earthquake  source.  Duration  of  ground  motion  and 
period  estimates,  based  on  design  earthquake  parameters  md  and  R^,  provide 
valuable  information.  If  a  design  response  spectrum  has  been  determined,  an 
artificial  time-history  record  can  be  generated  using  a  computer  program  devel¬ 
oped  by  Vanmarcke  at  M.I.T. 

Seismic  Hazard  Maps 

A  seismic  hazard  map  can  be  generated  for  a  region  by  dividing  the  region 
into  a  grid  of  points  and  repeating  the  analysis  treating  each  point  as  a  site. 
Contour  lines  are  then  drawn  for  various  levels  of  acceleration  corresponding 
to  a  certain  probability  of  exceedence  in  a  given  period  of  time.  A  hazard  map 
for  a  region  provides  Information  regarding  the  expected  hazard  at  different 
locations  in  the  region  as  well  as  valuable  input  to  the  formulation  of  seismic 
provisions  in  building  codes.  In  addition,  hazard  maps  can  be  used  in  the  seis¬ 
mic  risk  assessment  for  lifelines. 

Various  maps  for  the  United  States,  of  different  forms  and  containing  var¬ 
ious  modifications,  have  been  published  since  1947.  Although  certain  of  the 
early  maps  were  labeled  as  risk  maps,  they  actually  provided  no  estimate  of 
risk  or  frequency  of  occurrence  of  seismic  events.  This  Report  will  concern 
itself  with  only  the  more  recent  and  important  of  the  published  maps. 

Uniform  Building  Code 

Fig.  4.1  shows  the  zoning  map  adopted  by  the  Uniform  Building  Code,  In 
which  the  United  States  is  divided  into  five  zones  identified  as  zones  0,  1,2, 
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ZONE  1  -  MINOR  DAMAGE 
ZONE  2  -  MODERATE  DAMAGE 
ZONE  3  -  MAJOR  DAMAGE 
ZONE  4  -  GREAT  DAMAGE 


3,  and  4,  and  Indicates  the  different  expected  damages  for  each  zone.  This  map, 
although  widely  applied  In  engineering  practice,  contains  certain  shortcomings. 
The  map  does  not  take  Into  consideration  the  occurrence  frequency  of  earthquakes, 
and  therefore  cannot  make  any  distinction  between  regions  of  different  earth¬ 
quake  frequencies.  Thus,  an  assessment  of  seismic  risk  using  this  map  is  impos¬ 
sible. 

California  Division  of  Mines  and  Geology 

Greensfelder  (1973)  presented  a  hazard  map  for  California  in  which  contours 
of  maximum  credible  rock  acceleration  were  shown.  This  map  also  contains  certain 
shortcomings.  The  maximum  creditable  event,  a  figure  which  is  generally  ar¬ 
rived  at  based  on  geologic  evidence  using  a  geologic  timespan,  is  critical  for 
major  structures  where  the  consequences  of  failure  are  quite  serious.  However, 
in  regard  to  most  ordinary  structures  with  economic  life  spans  of  50  -  200  years, 
design  criteria  based  on  the  maximum  credible  event  would  be  over-conservative. 
Greensfelder' s  map  does  not  contain  any  information  regarding  the  probable  level 
of  loading  on  a  structure  or  the  frequency  of  occurrence  of  events. 

It  is  clear  that  in  order  to  assess  the  seismic  risk  associated  with  the 
hazard  level  suggested  by  a  map,  a  probability  level  must  be  associated  with  the 
hazard  parameter  Indicated  in  the  map. 

Milne  and  Davenport  (1969) 

The  first  map  for  the  United  States  which  Incorporated  probability  concepts 
was  published  by  Milne  and  Davenport  (1969).  The  two  maps  presented  by  them 
gave  the  return  period  In  years  for  acceleration  of  10%  gravity  and  the  accel¬ 
eration  in  percent  of  gravity  for  a  100-year  return  period. 
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Kiremidjian  and  Shah  0975) 

More  recently,  a  seismic  hazard  map  for  California  was  published  by 
Kiremidjian  and  Shah  (1975),  This  map  gives  peak  ground  acceleration  values 
corresponding  to  P[A  >  a]  =  0.1,  during  a  period  of  50  years.  The  return 
period  associated  with  this  level  of  probability  is  475  years. 

Alqermissen  and  Perkins  (1976) 

Probabilistic  estimates  of  maximum  acceleration  in  rock  in  the  contiguous 
United  States  were  presented  by  Algermissen  and  Perkins  (1976)  of  the  U.S.  Geo¬ 
logical  Survey.  The  map  shown  in  Fig.  4.2  indicates  the  relative  hazard  for 
various  parts  of  the  United  States.  The  contour  lines  for  horizontal  accelera¬ 
tion  in  percent  gravity  correspond  to  a  90%  probability  of  not  being  exceeded  in 
50  years.  The  seismic  hazard  analysis  conducted  to  arrive  at  these  estimates 
was  based  primarily  on  historic  data.  The  acceleration  attenuation  laws  employed 
in  the  analysis  for  the  Western  and  Eastern  United  States  were  shown  in  Fig.  2.14. 
Applied  Technology  Council  (1978) 

The  design  earthquake  ground  motion  and  elastic  spectrum  provisions  developed 
by  the  Applied  Technology  Council  (ATC),  associated  with  the  Structural  Engineers 
Association  of  California,  were  recently  publ ished  (5,  105).  The  two  parameters  em¬ 
ployed  to  describe  the  provisions  are  "effective  peak  acceleration"  (EPA),  and  "effect 
ive  peak  velocity"  (EPV).  These  parameters  are  used  as  scaling  factors  to  de¬ 
velop  a  design  response  spectrum  from  standardized  spectra.  The  parameters  do 
not  have  physical  meaning,  but  were  chosen  such  that  the  EPA  is  proportional  to 
a  spectral  ordinate  between  0.1-  to  0.5-second  periods,  and  the  EPV  is  proportional 
to  a  spectral  ordinate  of  approximately  1-second  period.  This  proportionality 
constant  is  2.5  for  a  damping  value  of  5%.  Thus,  from  a  design  spectrum,  it  is 

possible  to  calculate  EPA  *  S-,(0.1  to  0.5)  and  EPV  =  (1  sec)  #  jn  general, 

2.5  2.5 
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(FrS!  ATCigTS)”315  Prepared  by  Al9ermissen  a"«  Perkins  (1976) 


EPA  is  less  than  the  peak  ground  acceleration,  particularly  for  very  high  fre¬ 
quencies;  and  EPV  is  greater  than  peak  ground  velocity,  particularly  at  long  dis¬ 
tances.  The  duration  of  ground  motion  associated  with  the  EPA  and  EPV  values 
recommended  in  the  provision  is  between  20  and  32  seconds.  These  values  should 
be  decreased  for  shorter  durations  and  increased  for  longer  durations,  as  deemed 
appropriate. 

Seismic  hazard  analyses  were  conducted  using  the  Algermissen  and  Perkins 
(1976)  map  as  a  reference.  The  results  were  interpreted  as  a  map  for  EPA,  shown 
in  Figs.  4.3  and  4.4.  The  map  for  EPV,  shown  in  Figs.  4.5  and  4.6,  was  obtained 
by  introducing  modifications  to  the  map  for  EPA.  The  numbers  on  the  contour 
lines  in  Fig.  4.3  correspond  to  Afl,  the  value  of  EPA  divided  by  gravity.  The 
values  in  Fig.  4.5  correspond  to  Ay,  which  is  a  dimensionless  parameter  re¬ 
ferred  to  as  velocity-related  acceleration,  and  is  related  to  EPV  as  follows: 


EPV  (sec)  A>, 

12  0.4 

6  0.2 

3  0.1 

1.5  0.05 


The  contour  lines  correspond  to  a  90%  probability  of  not  exceeding,  in  a 
period  of  50  years,  the  level  of  risk  chosen  by  Algermissen  and  Perkins.  If 
the  EPA  chosen  for  design  differs  from  the  value  suggested  by  the  map,  the  pro¬ 
bability  of  not  exceeding  in  50  years,  corresponding  to  the  value  of  EPA,  can 
be  determined  from  Fig.  4.7. 
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ALASKA 


PUERTO  RICO 


FIG.  4.6  Contour  Map  for  Effective  Peak  Velocity-Related 
Acceleration  Coefficient  (From  ATC  1978) 


117 


ANNUAL 


PROBABILITY  THAT  MOTION  NOT  EXCEEDED  IN  50  YEARS 


Design  Spectrum 

Having  estimated  the  EPA  for  a  site,  the  design  spectrum  for  horizontal 
motion  can  be  determined  by  using  the  normalized  spectra  shown  in  Fig.  4.8 
for  three  soil  types.  Briefly,  refers  to  a  soil  profile  consisting  of  rock 
or  stiff  deposits  (less  than  200  ft*)  of  sand,  gravel  and  stiffer  clays;  S2  re¬ 
fers  to  deep  (greater  than  200-ft)  cohesionless  soils  or  stiff  clay  deposits; 
and  S3  refers  to  deposits  of  30  feet  or  more  of  clays,  sands  and  other  cohesion¬ 
less  materials. 

In  order  to  determine  the  design  elastic  spectrum  for  5%  damping,  the  appro¬ 
priate  spectrum  from  Fig.  4.8  is  multiplied  by  EPA  x  S,  where  S  =  1.0  for 
and  S,,  profiles,  and  S  =  0.8  for  S3  profiles.  The  design  spectrum  for  2%  damp¬ 
ing  can  be  determined  by  multiplying  the  design  spectrum  for  S%  damping  by  1.25. 
The  design  elastic  response  spectrum  for  vertical  motions  is  determined  by  multi¬ 
plying  the  horizontal  spectrum  by  0.67. 

The  procedure  described  is  applicable,  provided  that  Aa  and  Av  are  identical, 
as  in  the  region  of  the  innermost  contours.  Where  A  and  A  differ,  the  pro- 

d  V 

vision  states  that  "the  portion  of  the  response  spectrum  controlled  by  the  velo¬ 
city  should  be  increased  in  proportion  to  the  EPV  value  and  the  remainder  of  the 
response  spectrum  extended  to  maintain  the  same  overall  spectral  form." 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measurement  to 
metric  (SI)  units  is  presented  on  page  3. 
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SPECTRAL  acceleration 


FIG.  4.8  Normalized  Response  Spectra  Recomnended  for  Use 
in  Building  Code  (From  ATC  1978) 


PART  5:  CONCLUSIONS 


Seismic  hazard  analysis  provides  estimates  of  the  probability  that  a  de¬ 
sign  parameter,  P,  at  a  given  site,  will  exceed  the  value  p  within  a  certain 
period  of  time.  Essentially,  the  procedure  Involves  the  following  general 
steps. 

1)  A  compilation  of  geologic,  geophysical  and  selsmologlcal  data  regard¬ 
ing  the  seismicity  of  the  region  under  Investigation. 

2)  The  establishment  of  potential  earthquake  sources  within  the  region 
using  all  the  sources  of  information. 

3)  For  each  source,  the  determination  of  the  upper  and  lower  bound  magnitudes 
(intensities),  the  mean  annual  number  of  earthquake  occurrences,  and  the  rela¬ 
tive  frequency  distribution  for  earthquake  magnitude  (Intensity). 

4)  Selection  of  an  appropriate  attenuation  law  to  relate  the  seismic 
design  parameter  at  the  site  to  earthquake  magnitude  (Intensity),  and  distance 
to  source  of  energy  release. 

5)  Calculation  of  the  probability  of  exceedence  (or  non-exceedence)  of 
p  for  various  assumed  values  of  P  using  probabilistic  procedures. 

Calculation  of  the  input  parameters  required  for  the  analysis  involves 
numerous  uncertain! tes  which  can  be  accounted  for  and  incorporated  in  the  pro¬ 
bability  calculations.  It  must  be  noted  that  the  reliability  of  the  results  of 
seismic  hazard  analysis  depends  on  the  accuracy  of  these  Input  data.  Caution  Is 
also  advised  regarding  the  calculation  of  small  probabilities  which  are  associ¬ 
ated  with  the  tails  of  the  probability  distribution  for  the  parameter  P  used  in 
the  analysis  (38).  Proper  Interpretation  of  the  results  requires  the  applica¬ 
tion  of  professional  judgment  and  experience. 
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The  practical  application  of  seismic  hazard  analysis  In  most  cases  has  been  lim¬ 
ited  to  the  definition  of  a  design  seismic  parameter  corresponding  to  an  acceptable 
level  of  risk.  The  results  of  seismic  hazard  analysis,  however,  are  beneficial 
beyond  their  present  practical  use  when  applied  In  research.  Incorporating  seis¬ 
mic  hazard  analysis  results  into  an  overall  risk  analysis  can  identify  the  major 
seismic  parameters  contributing  to  the  risk.  Furthermore,  It  may  lead  to  alter¬ 
nate  and  more  realistic  procedures  for  analysis  In  earthquake  engineering. 
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APPENDIX  A:  SEISMIC  HAZARD  ANALYSIS  USING  PROBABILISTIC 
ATTENUATION  LAW 


McGuire  (1976)  presents  a  brief  description  of  the  theoretical  formulations 
used  in  his  computer  program  for  seismic  hazard  analysis  employing  a  probabilistic 
attenuation  law.  The  mathematical  expressions  presented  herein  are  taken  from 
McGuire  (1976). 

The  unconditional  probability  of  acceleration  A  exceeding  a  is  given  by: 

P[A  >  a]  =  /  J  P[A  >  a|m  and  R]  fM(m)  fR(r)  dmdr  (A.l) 

where  P[A  >  a|m  and  R]  is  the  complementary  cumulative  probability  of  a,  given 
an  earthquake  of  magnitude  m  and  distance  R;  fM(m)  and  fR(r)  represent  the 
probability  density  functions  of  m  and  R.  The  integration  is  carried  out  over 
the  entire  magnitude  range  of  engineering  interest,  and  over  the  various  sources. 
Assuming  an  attenuation  law  of  the  form 

‘  ■  C1  eC21"  (r  ♦  r/3  (A. 2) 

where  a  represents  the  mean  value  of  the  acceleration,  which  is  lognormally  dis¬ 
tributed  with  a  standard  deviation  <Jjna*  and  using  the  truncated  magnitude  re¬ 
currence  relationship  given  in  Eq.  2.10.  the  density  function  on  magnitude  is 
given  by: 

=  BKm1  exp(-B(m  -  m0)  )  mQ  <  m  <  m1  (A. 3) 


A1 


Using  Eqs.  A. 2  and  A. 3,  £q.  A.l  can  be  written  as: 


Is  evaluated  numerically  for  each  value  of  acceleration  a.  The  output  from 
McGuire's  program  includes  values  for  the  total  expected  number  of  events 
causing  intensities  or  accelerations  exceeding  some  specified  values,  as  well 
as  the  annual  probabilities  of  exceedence  associated  with  each  specified  event. 
The  probability  of  exceedence  during  a  time  interval  t  can  be  calculated 
from  the  computed  expected  number  of  occurrences  of  an  event  using  Eq.  3.5. 
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